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ABSTRACT 


This thesis is focused on solving one component of the proposed problem 
in the Global Trajectory Optimization Competition released by the Jet Propulsion 
Laboratory in late 2006. The goal is to find an optimal spacecraft trajectory to 
rendezvous with an asteroid in a group of asteroids. The analysis is conducted 
using a MATLAB application package for dynamic optimization called DIDO. In 
order to verify the selection results, one-to-one transfers between Earth and 
several asteroids are conducted. The selection process is applied to this group 
of asteroids. When the initial results do not meet the expectations based on the 
one-to-one transfers, a more thorough search for a global minimum is necessary. 
The gradual cost-constrained technique is used to progress from local minima 
toward the global minimum. The results are checked to satisfy the constraints as 
well as the necessary conditions for optimality. When the results are analyzed, 
feasible one-to-one rendezvous trajectories are found, however a sufficient 
selection process is lacking. There is a great deal of work remaining on this 


project, including the continued development of an asteroid selection procedure. 


THIS PAGE INTENTIONALLY LEFT BLANK 


vi 


TABLE OF CONTENTS 


I. INTRODUCTION viistactuanctiuindantdonedentaniudanideaadoatdeiedacadnnndaubdetedacuiuindenedebedentaviudand 1 
A. MOTIVATION 12 tivesdtcpriasstenupestiarneacttatsecendstinasecapevasdeepinanevapecssntastunctes® 1 

B. THE PROBLEM siscscisssscvcis saeccneccisectewnscassccctcasasols canveasceuesscseceusaceieotened 1 

C. HISTORICAL BACKGROUND ............ccsseeeeeeeeeeeeeeneeeeeeeeeeeeeeeeeneneeeeneees 3 

D. PSEUDOSPECTRAL METHOD .........::ccccecceeseeeeeeeeeeeeeeeeeeeeeeeeeeneeeeeeees 5 

E; GONGLUSION si isaiacstnstsaciiadaintaabotetcasteieatiwaateeescmolonnuieinaancoaneneintniennais 8 

ll. THE PROBLEM FORMULATION vais wtiscsitedinsainnsnnnsetetuecsonvtenauteswanseneninheaidutueans 9 
A. INTRODUGTION srincitecttiinninnnnhinaninna annem 9 

B. GOVERNING DYNAMIGS ivictiisicccidecces cacceincnattateacedssvnecasencnesideucecsareaaue 9 

C. OPTIMAL CONTROL PROBLEM ...........:ccs::eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 12 

D. SOLUTION METHOD ssestasiaiesstastoarensonsinnntaadiaveabainnscunatnenaaoninsaanntnas 13 

E. NECESSARY CONDITIONS FOR OPTIMALITY .........::::ccccccesseeeeee 15 

F. COMPUTER PROGRAMMING ....0....ccccccsseseeeeeeeeeeeeeeeeeeeneneeeeeeeeeeeeeeees 19 

1. Model Components for DIDO............:::::eeeeeeeeeeeeeeeeeeeeeeeeeeeeees 19 

2. SG ANA isda iiacecawecatecuutcasicacsatudeteccuu (any tuaslinscuuduny inavieausuaviiueadels 20 

G. FEASIBILITY-OF SOLUTION titvisnucnirsiennicvisesenianiiannestonweandationsinraewwtes 21 

ll. ONE-TO-ONE TRANSFERS satcsssetseatcetideniienineadivticoahivibientdebiteatieahinubiuaieaahts 25 
A. INTRODUCTION isssitcansatacaessecsan duntdenitradsintveddenivaitivucterddeuiendivusdsatease 25 

B. ONE-TO-ONE RENDEZVOUS RESULTS ........ccccccesseseeseeeeeeeeeeeeeeeeees 25 

C. VERIFICATION OF THE RESULTS) ...........::::eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 26 

1. Feasibility ANAlY SIS iiccnisastersaeciciiisteiadiabereddanioinasaraaneaeaceonnins 26 

2. Checking the Boundary Conditions.............::::cccccceesseseeeeeee 38 

D. CORRECTING THE MODEL .vi.cicccsscccsiecicaceeissescsseueenrssesieasssncieenonsies 43 

1. Problem Id@ntiticau On isiiicacicaiicacicacactcasacacteicacecteasacantcianciaaeina 43 

2. Problem Correction .........::::::ccccccceeseseeseeeeeeeeeeeeeeeeneeeeeeeeeeeeeees 43 

IV... ASTEROID. SELECTION s wstasccscuitscnceranencteess cansavenascteusscesitesseneaansnecsdeusrauinawane 55 
A. SELECTION LOGIG wisascedesactsedinssthdeseutletatashindbenstledinehietitisblaretahtutes 55 

B. SELECTION RESULTS sc caxcciyssntesaresiruanasuteinyeastaanusceledayeibisastedeioniaans 55 

C. GLOBAL OP TIMALUVY .ivisiciciccaccisticatcsveds sivenctenssiaacecdvaseventuanateudwaccdans 60 

1. DIDO and Global Minima .........cccccsssccseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 60 

2. Gradual Cost-Constrained Optimization Resulfts................ 61 

D. VALIDATION OF RESULT visaisnisiicsinsses ctiteniintinabnsesicutifotihesaduabewtteaahes 62 

V. CONCLUSIONS wasionitinascostcttesaatviatercteus cinvidtvastotesuadschivanitaldneeseuuensasteisteadians 67 
A. INTRODUCTION sisseccauectstcanccausunscctudeansiavsuansabedusntiansunssutuduandutuseavianss 67 

B. SUGGESTIONS ws ceicinsscesseeieindcineenssinmnenteinenertnnseusieunienntnendnnntenemnentsins 67 

C. FUTURE WORK oiotiseisctivnseastasescie deci teen teranenetaue tuwaduceneatentavaedecanenedeubncdens 68 

D. CONCLUDING REMARKG ...........:::::eeceeeeeeeeeeeeeeeeeeeeeeeeeessneseneeeeeeeeeees 69 
APPENDIX A: ASTEROID: DATA FILE saiisccicccsssciacsiccsticcssnstarvinaniensannesannsais 71 
APPENDIX B: M-FILES USED FOR FEASIBLE SOLUTION ............:::::eeeee 91 


vii 


LIST OF REFERENCES 
INITIAL DISTRIBUTION 


LIS E satsichiessccwisazinanbsupissudedptsanbsnsisiahtusbsiabiutiuaabcatinibeatinens 101 


viii 


Figure 1: 
Figure 2: 
Figure 3: 
Figure 4: 
Figure 5: 
Figure 6: 
Figure 7: 
Figure 8: 
Figure 9: 


Figure 10: 
Figure 11: 
Figure 12: 
Figure 13: 
Figure 14: 
Figure 15: 
Figure 16: 
Figure 17: 
Figure 18: 
Figure 19: 
Figure 20: 
Figure 21: 
Figure 22: 
Figure 23: 
Figure 24: 
Figure 25: 
Figure 26: 
Figure 27: 
Figure 28: 
Figure 29: 


Figure 30: 
Figure 31: 


LIST OF FIGURES 


Sample Feasible DIDO Solution ...................cccccseeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 22 
Sample Infeasible DIDO Solution ...................cccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 23 
Visual Feasibility Plot for Asteroid 906 with 15 Nodes...............:: 28 
Visual Feasibility Plot for Asteroid 907 with 15 Nodes..............::08 29 
Visual Feasibility Plot for Asteroid 908 with 15 Nodes.................:005 30 
Visual Feasibility Plot for Asteroid 909 with 15 Nodes.............:.:: 31 
Visual Feasibility Plot for Asteroid 910 with 15 Nodes..............:. 32 
Visual Feasibility Plot for Asteroid 906 with 60 Nodes..............::0 33 
Visual Feasibility Plot for Asteroid 907 with 60 Nodes..................:005 34 
Visual Feasibility Plot for Asteroid 908 with 60 Nodes...................000 35 
Visual Feasibility Plot for Asteroid 909 with 60 Nodes... 36 
Visual Feasibility Plot for Asteroid 910 with 60 Nodes............... 37 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

ASIGoid 906 Oris cea ce a eet ree ae aacaaeie tee canare teatecss 40 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

ASIGIOIO-OO7 OMDUS) c:cissden.ch cgecnpasteediocabeieatasedaneebs han oe Gamneebeee aaa 41 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

Asteroid: 90S Ons tacnd tenn redeemed ndaridndean en Bde Prien dad creak niads 41 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

ASIElOId DOO: OMNIS, osc caaccask cuss. ssepevagpeangies esas dengan bsbeh cst doaumehaassamcancb ons 42 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

ASteroid: 910 Olt rab tact ccc orctag lena nea ceedwaresteatinetlentatndies temworealegletsat dads 42 
Rescaled 30 Node Visual Feasibility Plot for Asteroid 906................. 46 
Rescaled 30 Node Visual Feasibility Plot for Asteroid 907 ................. 47 
Rescaled 30 Node Visual Feasibility Plot for Asteroid 908................. 48 
Rescaled 30 Node Visual Feasibility Plot for Asteroid 909................. 49 
Rescaled 30 Node Visual Feasibility Plot for Asteroid 910................. 50 
Rescaled 3-D Plot of Spacecraft Trajectory from Earth to Asteroid 

S| 8 | s ARR PRRNPE Rie ERP ER PARADE RAIPR Rn” SEMPRA? ESE Rr BST? ERAT MEPPEAAD: RAPE DERN? SADE RAT ORAPE RADAR ATER AES 51 
Rescaled 3-D Plot of Spacecraft Trajectory from Earth to Asteroid 

OF vente ence cate cseceeecvae cuseecosceus Cent cesterezecntcavecececeseceseeceeeusbcesecooece teense eeneiene 51 
Rescaled 3-D Plot of Spacecraft Trajectory from Earth to Asteroid 

DOB sco 2 disci ecas ne ene le cue deietledasdus.dss duadie dade es etuaduseiets 52 
Rescaled 3-D Plot of Spacecraft Trajectory from Earth to Asteroid 

12] OS Renner Ree ceoeR On eR en ty Et os EC OREC RRC RRC DER OORRN EEO MRE SEENDERIDER CORE SER CPER DEE yen CDEEC 52 
Rescaled 3-D Plot of Spacecraft Trajectory from Earth to Asteroid 

S18 iar recer em Panto eee ee a Vicar i ea eu ta Wr Pen ete Ne Ae Pa ol PORTE eee 53 
3-D Plot of Spacecraft Trajectory with respect to the Earth and 

Asteroid 907 based on the New Selection Technique .................:::005 57 
Visual Feasibility of Solution with Selection of Asteroid 907............... 58 
Feasibility of Selected Asteroid Utilizing New Selection Technique ... 59 
Demonstration of Global Minimum ...............::::ccccccceeeeeeeeeeeeeeeeeeeeeeeeees 60 


ix 


Figure 32: 
Figure 33: 


Visual Feasibility of Selected Trajectory Solution to Asteroid 909...... 
3-D Plot of Spacecraft Trajectory with respect to the Orbits of the 
Earth and Target Asteroid 909 0.0.0.0... eee eeeeeeee 


Table 1: 
Table 2: 


Table 3: 
Table 4: 
Table 5: 
Table 6: 
Table 7: 
Table 8: 
Table 9: 


Table 10: 


Table 11: 
Table 12: 


Table 13: 


LIST OF TABLES 


Spacecraft engine Specifications ................cccceccceeeeeeeeeeeeeeeeeeeeeeeeeeeeeeees 2 
Keplerian orbit elements of the Earth, J2000 heliocentric ecliptic 

POISON CS MANNE aerctcarnte ernie ela tec ede dtaradteineh essen tinsne teu ele atenutaunentaniindn! 2 
CONSLANtS: ANG. CONVESIONS: aia; caus cantvessceskeads cueteast coatseuaceskeauadesteouegeseentss 3 
One-to-One Transfer Costs and Visible Feasibility..................:::0 26 
Bootstrapped One-to-one Transfer Costs and Visible Feasibility ....... 27 
Boundary Conditions for the 60 Node Solutions. ................::::::::seeeee 39 
Rescaled One-to-One ReEndeZvous ..............:ccccccccceeeeeeeeeeeeeeeeeeeeeeeeeees 44 
Boundary Conditions for the Rescaled One-to-One Rendezvous ...... 45 
DIDO Output for a Model with Five Asteroids .................ccccceseeeeeeeeeeees 55 
DIDO Output for a Model with Five Asteroids with new Selection 

ESCM setendSest ible ccdcecceese detcdee ob dcwad dubcowidedbcceddevhe cuadenblouadecbeedce Aaceent 56 
Endpoint Constraints for new Selection Technique ................:::::cccee 57 
Selection Among Five Asteroids (906-910) with the Gradual Cost- 

Constrained : FEGMIMIQU CS 2s25uc5cs-hiiec case euencasbaesrahennaseaaahinas exe rbieetereuenacsh ene 62 
Boundary Conditions for the Selected Trajectory Solution to 

Fats] (<1 f6)[0 Wo) 01° tant ane ean eae eee Ao ee ee ee ee eee ee eee ee ee 63 


xi 


THIS PAGE INTENTIONALLY LEFT BLANK 


xii 


ACKNOWLEDGMENTS 


It would be inappropriate if | did not acknowledge, first and foremost, the 
assistance and encouragement that | received from Dr. Pooya Sekhavat and 
Professor |. M. Ross. The GNC Lab has been a great home for me over the past 
several months. | greatly appreciate all your patience and support while 
conducting this research. | would have never been able to complete this work if it 


weren't for your continuous sharing of insight and assistance. 


| also must thank my beautiful fiancée, Alison. Your love and support of 
me, even when | doubted myself, was undying and steadying. | appreciate all 
that you are to me and thank you for your undying patience while completing this 


work. 


Thanks must also go out to all of the Ensigns who have supported me in 
so many ways during our year together in Monterey. My time here would not 


have been the same if it weren’t for all of you. 


xiii 


THIS PAGE INTENTIONALLY LEFT BLANK 


Xiv 


I. INTRODUCTION 


A. MOTIVATION 


On 5 October 2006, the Outer Planets Mission Analysis Group of the Jet 
Propulsion Laboratory (JPL) released the second problem in the Global 
Trajectory Optimization Competition (GTOC2). The problem specified is a low- 
thrust trajectory design for a multiple asteroid rendezvous. The goal is to 
maximize the ratio of the final mass of the spacecraft to flight time. There are 
four clusters of asteroids to be met in sequence. The spacecraft is to 
rendezvous with one asteroid in each of four asteroid clusters. The order in 
which the spacecraft visits the clusters is inconsequential. This problem was 
chosen by JPL based on four criteria: 1) global optimization consists of several 
local optima over a large launch window, 2) unusual objective function or 
constraints, 3) the problem can be solved in a reasonable amount of time, and 4) 
solutions to the problem can be easily verified.1_ The goal of this research is to 
solve part of this problem, the selection of an asteroid in a group of asteroids 


which provides the optimal trajectory for rendezvous. 
B. THE PROBLEM 


This thesis focuses on the selection of an individual asteroid in a cluster 
with which to rendezvous, such that the cost is minimized; this work, therefore, is 
only addressing one component of the GTOC2 problem. The problem is 
simplified by limiting the forces acting on the spacecraft to only the gravitational 
force of the sun and the thrust provided by the spacecraft’s engine. As such, 
there can be no gravity assists used in the solution of the problem. The 
propulsion for the spacecraft is provided by a thruster which can be turned on 


1 Anastassios E. Petropoulos, “Problem Description for the 2"° Global Trajectory Optimization 
Competition,” Outer Planets Mission Analysis Group, Jet Propulsion Laboratory, 6 November 
2006. 
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and off at will. The direction of the propulsion is unconstrained. The 


specifications for the engine are provided in Table 1.2 








Table 1: Spacecraft engine specifications 
Constant Specific Impulse 4000 s 
Maximum Thrust Level 0.1N 














There is a twenty year window for the launch date, which must fall between 2015 
and 2035. The spacecraft is to launch from the Earth with a hyperbolic escape 


velocity of 3.5 km/s. The direction of this velocity is unconstrained. 


According to the problem statement, the Earth and asteroids can be 
assumed to follow Keplerian (conic) orbits about the Sun, eliminating 
perturbations and thus further simplifying the problem. The Keplerian orbital 


elements for the Earth are provided in Table 2.3 




















Table 2: Keplerian orbit elements of the Earth, J2000 heliocentric 
ecliptic reference frame 
Semimajor axis, a (AU) 0.999988049532578 
Eccentricity, e 0.01671681 163160 
Inclination , i (deg) 0.0008854353079654 
Longitude of ascending node, 0 (deg) | 175.40647696473 
Argument of perigee, w (deg) 287.61577546182 
Mean anomaly at epoch, Mo (deg) 257.60683707535 
Epoch, to (MJD) 54000 








2 Anastassios E. Petropoulos, “Problem Description for the 2"" Global Trajectory Optimization 
Competition.” 


3 Ibid. 





The constants and conversions which are to be used in solving this problem are 


provided in Table 3.4 


























Table 3: Constants and conversions 

Gravitational parameter of the Sun, yp (km*/s*) | 1.32712440018 x 10" 
Astronomical unit (AU), (km) 1.49597870681 x 10° 
Standard acceleration due to gravity, g (m/s?) 9.80665 

Day, (s) 86400 

Year, (days) 365.25 

00:00 01 January 2015, (MJD) 57023.5 

24:00 31 December 2035, (MJD) 64693.5 





The data for the asteroids are provided in an ASCII text file, containing the 
Keplerian orbital elements (see Appendix A). The following data are specified for 
each asteroid: 1) a unique ID number of the asteroid; 2) semimajor axis in AU; 3) 
eccentricity; 4) inclination in degrees; 5) longitude of the ascending node in 
degrees; 6) argument of periapsis in degrees; 7) mean anomaly in degrees; 8) 
epoch, in MJD, at which mean anomaly is given; and 9) the group number of the 
asteroid. Just as they are for Earth, the orbital elements of the asteroids are 
given in the J2000 heliocentric ecliptic frame.® From all of these provided values, 
constants and constraints, the dynamics of the problem can be developed. This 
will be shown in Chapter II. 


C. HISTORICAL BACKGROUND 


The motion of celestial bodies has been studied since the onset of 
civilization. In many ways, the level of development of a civilization can be 


determined by their knowledge of the workings of the heavens. The ability of 





4 Anastassios E. Petropoulos, “Problem Description for the 2"° Global Trajectory Optimization 
Competition.” 





humans to predict the motion of objects in space has become quite sophisticated 
through time, taking into account minuscule forces which affect the objects in 
minute, but distinguishable ways. For this research, the lack of perturbations to 
the orbits permits the simple use of the equations generated by Johannes Kepler 
and Isaac Newton to predict the motion of the bodies. These two 
mathematicians generated equations that could predict the position of celestial 
bodies based on careful observations. Kepler used the data collected throughout 
Tycho Brahe’s life to come up with his three laws; 1) the orbits of the planets are 
ellipses with the Sun at one focus, 2) the line joining a planet to the Sun sweeps 
out equal areas in equal times and 3) the square of the orbital period is directly 
proportional to the cube of the mean distance between the Sun and the planet.6 
A century later, Newton would use Galileo’s studies of falling objects to create his 
universal law of gravitation. This law states that the force of gravity between two 
bodies is directly proportional to the product of their masses divided by the 
square of the distance separating them. These simple ideas, developed in the 
16" and 17" centuries, will be used to develop the dynamics of the problem at 
hand. 


The desire for optimized solutions traces its roots back to the 17™ century 
when scientists and mathematicians such as Galileo Galilei and Pierre de Fermat 
posed problems which necessitated optimization. In order to solve these 
problems, Isaac Newton invented calculus of variations.’ Johannes Bernoulli 
was also instrumental in the initial surge for optimization as he used a discrete- 
step method to solve the brachistochrone problem in 1697. Many other great 
minds, such as Leonard Euler, Jean Louis Lagrange, Adrien-Marie Legendre 


(who first introduced multiple variable variations), Karl Gustav Jacob Jacobi and 





5 Anastassios E. Petropoulos, “Problem Description for the 2" Global Trajectory Optimization 
Competition.” 


6 Jerry Jon Sellers, Understanding Space: An Introduction to Astronautics (Boston: McGraw 
Hill Custom Publishing, 2004), 36-39. 


7 Arthur E. Bryson, Dynamic Optimization (Menlo Park, CA: Addison-Wesley Longman, Inc., 
1999), 403-405. 
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William Rowan Hamilton, made contributions to the field in the 18'" and 19" 
centuries. Numerous others have continued the rapid development and 
advancement of optimal control theory up to present day including Lev 
Pontryagin and Richard Bellman. 


D. PSEUDOSPECTRAL METHOD 


One of the most modern control theories is based on pseudospectral 
approximation. The pseudospectral (PS) approach can be implemented to solve 
dynamic optimization problems. It is used to solve the optimal control problem in 


this work. Here, we take a closer look at the theory behind it. 


Several pseudospectral methods can be used to solve optimal control 
problems. In this work, the Legendre Pseudospectral Method is utilized to 
numerically solve the dynamic optimal control problem. The advantage of the 
Legendre PS method is the ability to vary the spacing between data points.2 By 
placing the nodes at the roots of the Legendre polynomials, the integration of any 
polynomial of degree 2N-1 can be exactly evaluated. The Legendre PS method 
places the nodes at the zeros of the derivatives of the N-1 degree Legendre 


polynomials P,_,(t). The values at the endpoints are reserved for i=1 and i=N. 


The derivation of matrix D‘ necessary for implementation of the Legendre PS 


method is now possible through differentiation of the interpolation polynomial.1° 


XH) = X(t b(t) (1.1) 
X(t) = Dt), (t,) = DD); (1.2) 
D; = $(t,) (1.3) 





8 Arthur E. Bryson, Dynamic Optimization, 403-405. 


9 |. M. Ross and F. Fahroo, “Pseudospectral Methods for Optimal Motion Planning of 
Differentially Flat Systems,” (paper presented at IEEE Conference on Transactions on Automatic 
Control, Las Vegas, NV, Dec. 2002). 


10 |. M. Ross and F. Fahroo, “Legendre Pseudospectral Approximations of Optimal Control 
Problems,” Lecture Notes in Control and Information Sciences (2003). 
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The next step requires the interpolating polynomial to be in terms of the 
Legendre polynomial. The intermediate function is defined as: 


N 
z(t)=[[t-t, (1.4) 
i=1 
The derivative evaluated at a root, j, results in the following relationship. 


At) =(t,-t,)...(t) - t4)(t) -— tha).--(t) - ty) = [It -t, (1.5) 


Now, the Lagrange interpolating function is expressed in terms of the 
intermediate function (Eq. (1.4)).11 


é(t)=[] ie = Al) (1.6) 





The locations of the nodes {t,,...,fy_,} are chosen to be the roots of the derivative 
of the Legendre polynomials (Pat) It can be seen that the intermediate 
function has the same roots, with additional roots at the endpoints (1 and -1).12 

Z(t) =(t-£,)P'y s(t)(t= ty) = (t+ 1)P 'y (t)(t-1) = (8 -1)P y(t) (1.7) 
The Legendre differential equation is shown below. 


ia 
at 
Applying the Legendre differential equation to the intermediate equation yields 


[ (t? -1)P'y_4(t) ] = M(N-1)P, 4(t) (1.8) 


the following equation. 
2(t) = N(N-1)Py,,(t) (1.9) 
Combining Eqs. (1.6), (1.7), and (1.9) results in the following. 


(t? -1)P'y_,(t) 


1) = ENN 1)P,.AE) (1.10) 


11 |. M. Ross and F. Fahroo, “Legendre Pseudospectral Approximations of Optimal Control 
Problems.” 


12 |. M. Ross and F. Fahroo, “Pseudospectral Methods for Optimal Motion Planning of 
Differentially Flat Systems.” 
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Evaluating the derivative of (1.10) at t= t, will results in the following equation.'% 


1 2a PD alt ear 





0(t,) = N(N-1)P\,(t)| tt, i,t, (i, -t,)° 


Now, the Legendre differential equation (Eq. (1.8)) is simplified, 


N(N —1)P,, 4(t) OT (2 —1)P y(t) ] = 2tP ty (0) + (PP "yd (1.12) 


~ at 
And combined with Eq. (1.11) as follows. 14 


1 N(N-1)Py_(t)) _ (t? -1)P'y_(t) (1.13) 
(t;) i-f, (f, -t | 





Since (t? —1)P',_,(t,)=0, the differentiation matrix, D, 


;, can now be produced as: 





=4 = Py_s(t;) 
i OM Bb) ee 


The values which fill the differentiation matrix require further manipulation and 
analysis. For example, L’Hdpital’s rule must be applied to Eq. (1.13) so that the 


t;-t, terms can be eliminated from the denominator. From Eq. (1.14), the 


components of D; that form the N x N differentiation matrix DN, are:15 








P(t) 1 j 
Plt) Gat, 
N(N—-1) : 
Da ge ee! (1.15) 
N(N-1) adeG: 
4 
0 otherwise 


This completes the brief explanation of the numerical technique used in 
the Legendre Pseudospectral Method. A complete understanding of this theory 


13 |. M. Ross and F. Fahroo, “Legendre Pseudospectral Approximations of Optimal Control 
Problems.” 


14 Ibid. 


15 |. M. Ross and F. Fahroo, “Pseudospectral Methods for Optimal Motion Planning of 
Differentially Flat Systems.” 
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is not necessary for this work, as these calculations will be completed by DIDO, a 
MATLAB based application package authored at the Naval Postgraduate School. 


This will be explained in greater depth in the subsequent chapter. 
E; CONCLUSION 


The problem to be solved in this research is straightforward. The 
equations of motion will be derived based on simple two body orbital mechanics. 
There are several assumptions that are made which simplify the solution process 
(e.g. no gravity assists). The solution, however, is in no way simple to obtain. 
Optimal control theory will be applied to the dynamics governing the motion of 
the bodies such that the cost of the mission can be minimized. This problem is 
unique in that the minimized cost is not just the mass of the propellant or the time 
to travel the trajectory. Instead, the goal is to maximize the ratio of final mass 


and travel time, implementing both mass and time optimization. 


ll. © THE PROBLEM FORMULATION 


A. INTRODUCTION 


When JPL posed GTOC2, it was intended to be a problem that could be 
solved in a reasonable period of time by a team of experienced researchers. 
Therefore, the astrodynamics were simplified to a two body problem with only the 
gravity of the sun acting on the bodies. 


B. GOVERNING DYNAMICS 


The equations of motion of the problem stem from Newton’s laws. The 
only force acting on the Earth and the asteroids is the gravitational pull of the 
sun. Therefore, by combining Newton’s second law, 2F =™m.-a, with his law of 


-G-M-m_—-yu-m 
re fe 





gravitation, F= the equations of motion of the Earth and 


asteroids can be determined:16 


a ae =0 
J+ wz =0 (2.1) 
2+ =0 

r 


These are the two-body equations of motion in Cartesian coordinates (x,y,Z), 


where r=./x* + y*+z°. The gravitational constant, wv , in this case is that of the 


Sun, the value of which can be found in Table 3. 





16 Daryl G. Boden, “Introduction to Astrodynamics,” in Spacecraft Mission Analysis and 
Design, ed. Wiley J Larson and James R. Wertz (El Segundo, CA: Microcosm Press, 1999), 132. 
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The solution to the generic two body equation (= =0) is known to be 


oo 


the polar equation of a conic section.1” 


a-(1-e”’) 


~ (1+ e-cos(v)) SA) 





The velocity of an object in an elliptical orbit is determined by the following 
equation.18 

a sa (2.3) 

roa 

At this point, the complexity of converting between the elliptical coordinates (in 
which the orbital elements are presented in) and the Cartesian coordinates (in 
which the equation of motion are derived), begins to become apparent. In order 
to convert between these reference frames, the following relationships are 


used. 19 


r{[cos(6 + @)-cos(Q) - sin(@ + @) -cos(/)-sin(Q 


) 
r[cos( + @)- sin(Q) + sin(@ + @) -cos(/)-cos(Q)] 


r[sin(@ + @)-sin(i)] 
a 
[- 


[cos(@ + @—7)-sin(i)] 


(2.4) 


v|—sin(6 + @— v)-cos(Q)-cos(6 + @— v)-cos(/)-sin(Q 


x 
y 
Zz 
Vy ) 

Vy sin(9 + @— y)-sin(Q) + cos(6 + w— v)- cos(/)- cos(Q)| 
V, 


=V 
V 


These equations requite the flight path angle (7). This angle can be found from 


the following equation: 


e-sin(@) 


1+ e-cos(6) ee) 


tan(y) = 


which requires the value of the true anomaly (@). True anomaly can be found 


from: 


17 Jerry Sellers, Understanding Space, 136-137. 
18 Ibid. 


19 Anastassios E. Petropoulos, “Problem Description for the 2™ Global Trajectory 
Optimization Competition.” 
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tan( 5 = 8 tan{ 2) (2.6) 
2 1+e 2 


This equation contains the eccentric anomaly, E, whose value is found by using 
Kepler’s Equation, shown below: 

M=E-e-sin(E) (2.7) 
Kepler’s equation is a nonlinear equation that must be solved in an iterative 
process to calculate eccentric anomaly.29 It also requires the knowledge of 
mean anomaly, M, calculated from the following equation: 


mM, = | (4) (2.8) 


Where Mb is the initial mean anomaly of the object, to is the initial time when the 
object is at Mo, a is the semimajor axis of the orbit of the body and yw is that of 
the Sun. Eqs. (2.2) through (2.8) yield the Cartesian coordinates of the bodies in 


motion. 


The only difference between the equations of motion for the spacecraft 
and those of the celestial bodies is the addition of the spacecraft thrust force, T . 


This force can be projected into its Cartesian components. 


T=J/724+7,7 +72 <0.1N (2.9) 


The spacecraft equations of motion include this additional force as follows:21 


x x 

r m 

- 
jtpyaes (2.10) 

r m 

Zs 

r m 


The mass of the spacecraft will decrease as fuel is used to create thrust. 
Therefore, the last equation of motion is that of the mass flow rate of propellant: 


20 Anastassios E. Petropoulos, “Problem Description for the 2"° Global Trajectory 
Optimization Competition.” 


21 Ibid. 
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m= ae (2.11) 

Having developed the relationships which govern the motion of the Earth, 
asteroids and the spacecraft, the dynamics of the problem is established. These 
dynamic equations will constrain the solution of the optimal control problem and 


will be included in the optimal control problem definition. 
C. OPTIMAL CONTROL PROBLEM 


In this section the optimal control problem is formulated. We will establish 
the notation and methodology that will be used to solve the problem. The control 


function, u*(t), and the state function which results from the controls, x* (tf), 


must be established such that the Bolza cost functional is minimized. The Bolza 


cost functional is shown below.22 
J(X(-),U(-), tf) = E(x(t,),t,) + [ Flx(t),u(t), tat (2.12) 


In the Bolza cost functional, x<¢R” and uc R” are constrained by the 
differential constraint x = f(x(t),u(t)) where te<[t,,t,]. The boundary conditions 
can then be specified as:28 


@)(X(t,),f))=0 @ <€R?° and p<n+1 (2.13) 
e(x(t,),t,)=0  e ¢R? and q<n+1 


For this problem, the Bolza cost functional to be maximized is the ratio of final 
mass to flight time. 


J, = (2.14) 





For this research, the cost function will be the inverse of this value. This will 
allow for minimization of the cost function. Therefore, the Bolza cost functional 
which will be used for this analysis is: 





22 Arthur E. Bryson, Dynamic Optimization. 
23 Ibid. 


AZ 


jg (2.15) 





Based on the dynamics derived in Section B of this chapter, the state and control 
variables are: 

xeR’ {%Y.Z,M,V,,V,,Ve Bae 
ueR?’ (Tahal) DT 4+T/ +T al; VP +y+Z7 = rt , 


With knowledge of the governing dynamics of the problem (developed in Section 
B), and the cost function (Eq. (2.15)), the problem formulation for the spacecraft 


rendezvous is developed as follows: 





Minimize J[x(-),u(-),t; ] = aes 
Mm, 
Subject to X=V, 
y=V, 
Z2=V; 
. 7 
m=—— 
gl. 
Pode 
x m are 
hy y 
Mom AB 
ents Z 
Vv, =—— Ha 
mor 


(x6: Vorzo: MgsVy Vy Ve, ) = (Xess ihe dese rl 500 KQsV earth,» Vyearth,? V earth, ) 
v v )(2.17) 


(x ft Yr : Z; ! Vy, Vy, 4 Vv, ) = ( x, asteroid; ? y asteroid; ? Z asteroid, 1 V vasteroid, 1” yasteroid;? “ zasteroid; 


This is the problem that must be formulated and fed into DIDO to generate an 


optimal trajectory. 


D. SOLUTION METHOD 


Optimal control solution methods have been divided into two categories: 
indirect and direct methods. In indirect methods, Pontryagin’s Minimum Principle 
13 


is used to derive the necessary conditions and to obtain the optimal trajectory. 
Indirect methods are known to require rigorous efforts. Direct methods discretize 
the optimal control problem into a parameter optimization problem. The result is 
a nonlinear programming problem which can be solved based on known 
nonlinear programming solution methods. The MATLAB-based software 
package, DIDO, will be used to solve this problem. DIDO uses a direct method 
to solve optimal control problems. 


DIDO is a reusable software package which utilizes the Legendre 
Pseudospectral method to numerically solve partial differential equations. The 
Legendre Pseudospectral method is also commonly used in fluid mechanics. 
The uniqueness of the method comes from its use of global orthogonal 
polynomials as trial functions.24 It requires that the time domain be discretized 
into a set of Legendre-Gauss-Lobatto (LGL) points which will henceforth be 
referred to as nodes. The values of the state and control variables at the nodes 
are the unknowns while the Lagrange polynomials act as the trial functions.25 
Based on the Lagrange polynomials, the nonlinear state equations are 
transformed into nonlinear algebraic functions. The result is a nonlinear 
programming problem which can be solved using a sparse numerical optimizer. 
The relationship between the Lagrange variables and the costate variables 
allows the costates to be determined at the nodes based on the Covector 
Mapping Theorem.26 The Covector Mapping Theorem suggests that given any 
feasible solution to a problem, there exists a feasible solution to the augmented 





24 |. M. Ross and F. Fahroo, “User’s Manual for DIDO 2001: A MATLAB Application Package 
for Dynamic Optimization” Tec. Rep. NPS Technical Report AA-01-003, Department of 
Aeronautics and Astronautics, (Naval Postgraduate School, Monterey, CA, 2001). 


25 Ibid. 


26 Qi Gong, |. M. Ross, Wei Kang and F. Fahroo, “On the Pseudospectral Covector Mapping 
Theorem for Nonlinear Optimal Control” (paper presented at the 45" IEEE Conference on 
Decision & Control, San Diego, CA, Dec 13-15, 2006). 
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Karush-Kuhn-Tucker conditions, within a feasibility tolerance.2”? The Karush- 


Kuhn-Tucker conditions state the following.28 


<0 h,(u,t) = Ar(t) 
=0 if hr(t)< h(u,t)< Ae (t 
i ; (t) < A,(u,t) < h; (t) (2.18) 
>0 h(u,t)= A(t) 
unrestricted hr(t) = h?(t) 


where yz, is the covector associated with path constraints and h> and h” are the 
lower and upper bounds of the path constraint function (h,). Knowing the values 


of the costates allows for the validation of the optimality of the provided solution, 
as explained in the next section.22 More information on Pseudospectral methods 
can be found elsewhere. 


E. NECESSARY CONDITIONS FOR OPTIMALITY 


In this section, the the necessary conditions for optimality will be 
developed. These conditions are used to verify the optimality of a solution. The 
endpoint Lagrangian is known to be: 

F=Eiv'e (2.19) 
Where E is the endpoint cost function, v is the covector associated with endpoint 
constraints and e is the endpoint constraint function. The endpoint constraints 
present in the problem are the initial position, velocity and mass of the Earth and 


the final position and velocity of the asteroid. 


=| XypVos!Mys X 7p Mean» Mass, | (2.20) 





27 Qi Gong, I. M. Ross, Wei Kang and F. Fahroo, “On the Pseudospectral Covector Mapping 
Theorem for Nonlinear Optimal Control.” 


28 Qi Gong, |. M. Ross and Wei Kang, “A Pseudospectral Method for the Optimal Control of 
Constrained Feedback Linearizable Systems” (paper presented at the 2005 IEEE Conference on 
Control Applications, Toronto, Canada, August 28-31, 2005). 


29 Ibid. 
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The Hamiltonian is defined by the following relationship: 
H(A, x,u,t) = F(x,u,t) + 2"F(x,u,t) (2.21) 
where t—>2e¢R*“is a Lagrange multiplier function, also known as a costate.3° 


The costate must satisfy the differential equation known as the adjoint equation: 
A =-— (2.22) 
The state dynamics joined with the adjoint equation form a Hamiltonian system. 
For a control function u(-) to be optimal, it is necessary that u globally minimize 
the Hamiltonian for any te|t,,t,|.21 Ultimately, a candidate solution for optimal 
control results from the minimization of the Hamiltonian with respect to u if 2 and 


X are held constant. This is known as the Hamiltonian Minimization Condition. A 


candidate control function, u*(-), can be generated at each t. This control 
solution is Known as the extremal control. 
(A,x,t) Ke u* (2.23) 
Pontryagin’s Minimum principle develops the necessary conditions for 
optimality of the control u* based on the following equations.2 


Hamiltonian Minimization > H(A*, x*,u*,t) < H(A*, x*,u, t) (2.24) 
Adjoint equations > i* = = (2.25) 
x 

ry rT 0G, 
Initial transversality > A(t,) =—-v ae (2.26) 

Xo 
Terminal transversality > A(t,) = cls +v' oes (2.27) 

OX, OX, 

Hamiltonian Value —> H[t,] + & tyr < =0 (2.28) 


Based on the problem formulation in Eq. (2.17) and F(x,u,t)=0 the Hamiltonian 


then becomes: 





30 |. M. Ross, Control and Optimization: An introduction to Principles and Applications, 
Electronic Edition, Naval Postgraduates School, Monterey, CA, December, 2005. 


31 Ibid. 
32 Arthur E. Bryson, Dynamic Optimization. 
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= ie 
H(A, X,u,t) =A,V, + AV, + AV, + An p +2, (4 Lu ht Cae hla Ca 
x y m ze m 


g:l eo 


ais m r 


(2.29) 

Although it does not directly appear in Eq. (2.29), all of the terms are time 

dependant. This is a result of the time dependence of mean anomaly (Eq. (2.8)) 
as developed as a part of the equations of motion derived in Section B. 

oH 

ot 

The end-point Lagrangian can be defined, allowing the final value of the 


+0 (2.30) 


Hamiltonian to be determined. This information is useful when determining the 
validity of a solution produced by DIDO. 
oe 


Alles, =° (2.31) 


Since the endpoint conditions in this problem are also time dependant through 
the propagation of mean anomaly, E is dependent on time: 


GE _ OE " a(v'e) 1° Av'e) 








aes (2.32) 
at, at, at, m, dt 
T 
Gigjessa 2 (2.33) 
m, ot, 


The adjoint equations can also be developed from the Hamiltonian (see 
Eq. (2.25)). The adjoint equation defined earlier is applied to the Hamiltonian 


Lagrangian, H, defined as: 


A(z,4,X,0,t) = H(A,x,0,t)+ 7” -h(a) (2.34) 
where: 
__ [h(u)] |72+7?+72-T? =0 
h(u) = Sal , (2.35) 
AG)). ery ez =r =0 


The adjoint equations are therefore: 
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pageants 
Ox 
OH 
i a ey 2 [lyX 
. OH a 
Ay = dy 7 A,, 3 2LyY 
oH a 
a ae ey We ea 702 
Z 6z Vv, r° Hy 
: T 
See : Meg: —S +A, , iB 
om x om yom zm 
Piatra, 
* Ov, 
2 on 
Ov, 
= OTe 5, (2.36) 
7 OV 


Using Eq. (2.36), the dual variable outputs can be analyzed and used for further 
verification of the solution. 
Another option for analytically evaluating the optimality of a solution is 


through the analysis of the transversality conditions (Eqs. (2.26) and (2.27)). The 
endpoint constraint function of this problem is affected by changes in the states. 








A(t) = vi Le V, es V> 0g, Ve 085 (2.37) 
OX, OX, OX OX, 
A(t) = Oe 4 v7 2 = 04 y, 24 y, OE 4 y, & (2.38) 
OX, OX, OX, OX, OX, 


Solving these equations requires knowledge of the values of the covector. 


Unfortunately, these values are unknowns. 
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F. COMPUTER PROGRAMMING 


This section is intended to familiarize the reader with the mathematical 
model developed to solve of this problem. In order to utilize DIDO, the problem 
formulation must first be fed into it in the proper way. If the problem is not 
formulated and coded properly, the solution will not be correct as DIDO will be 
solving a different problem altogether. A great deal of the focus of this work has 
gone into the development and validation of the model. 


1. Model Components for DIDO 


In order to properly code the problem for DIDO, five MATLAB m-script files 
are created to define various components of the problem formulation. First, the 
cost function is defined in a cost m-file. This file contains the definition of the 
cost function which was developed in Eq. (2.15). Second, the dynamics of the 
problem is defined. This is done in an m-file which contains all of the dynamic 
equations developed in Section B. A third file needs to be created which states 
the end-point constraints based on the problem boundary conditions. This file is 
referred to as the events file. The events ensure that the spacecraft starts at 
Earth with the desired escape velocity and ends at an asteroid within the desired 
rendezvous tolerances. The asteroid selection logic is also included in this file. 
The fourth file is the path file and it is where the constraints along the path of the 
spacecraft are enforced. The path constraints are defined at the nodes based on 
the primal states and controls.%3_ The fifth file is the main file that calls the other 
four. This main file defines any constants and scaling factors and calls DIDO. 
The main file is necessary for the generation of all the necessary inputs for 
DIDO, such as the specifications of the nodes, bounds on the variables and the 
initial guesses. DIDO then takes all of the provided data and runs a search for 
an optimal control solution. When the model is properly formulated, and the 





33 |. M. Ross and F. Fahroo, “User’s Manual for DIDO 2001.” 
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problem statement accurately delineated, DIDO will provide the optimal solution. 
All of the files mentioned above are provided in Appendix B. 


2. Scaling 


It has already been stated that this is a low-thrust optimization problem, as 
such, the thrust force is quite small. The distances in question, on the other 
hand, are quite large. The problem is being solved numerically through the use 
of an embedded iterative optimization program. This introduces very important 
factors for consideration, round-off and truncation errors. Due to the large 
differences in the magnitudes of the values being used in this problem, it is 
necessary to scale all of the variables so that they will have equal importance in 
the analysis. The use of canonical units can be implemented to normalize the 
values of the variables and greatly diminish the effects of round-off error on the 
final solution. A common “Distance-Unit” (DU), “Time-Unit” (TU), “Mass-Unit” 
(MU), “Velocity-Unit” (VU), and “Force-Unit” (FU) can be generated based on the 
scale of the problem at hand.°4 This form of units is extremely valuable in 
solving the problem using a numeric approach. The variables must, however, be 
related through given values. An example of the development of a set of 


canonical units is shown below. 


Let DU=10-AU and Let MU=1000 kg 








DU* 1 rd 
T= based on Kepler's 3" law (2.39) 
L 
ruUS= es based on Newton's Law of Graviation (2.40) 


These units can be developed further into other values such as the velocity-unit 
(VU=DU/TU). It is important to balance all equations so that all of the variables 


used are equally scaled. Great difficulty can come from poorly scaled or 





34 Jerry Jon Sellers, Understanding Space, 699-700. 
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balanced equations.%5 If proper care and attention is paid to the scaling of a 
problem these difficulties can be avoided and the problem will be significantly 
simplified. 


G. FEASIBILITY OF SOLUTION 


The solution found by DIDO should be checked for both optimality and 
feasibility. In order to verify the feasibility of the solution, the control trajectory 
produced by DIDO should be propagated through the dynamics of the problem. 
The propagated states at the nodes should match the values generated by 
DIDO. These values can be plotted together to visually check the feasibility of 
the solution provided by DIDO, i.e., if the propagated solution passes through the 
nodes the solution is visually feasible. DIDO calls for the specification of the 
number of nodes that the solution will be based on. As a first guess, to speed 
the solution process, 15 nodes were used to solve the problem. This result was 
then inserted into DIDO (bootstrapped) as a guess for a solution with 30 or more 
nodes. This higher number of nodes increases the accuracy of the result 
allowing for better analysis of the results. In order to verify that the plot of the 
states was generated based on the proper propagated control solution, the 
interpolated control variables that are used for state propagation are checked 
against their node point values from DIDO as well. Figure 1 and Figure 2 
illustrate this consideration. It is important to note that a feasible solution in no 
way means an optimal solution has been generated. Rather, it only means that 
the state/control solution generated by DIDO satisfies the dynamics of the 
problem. 


To conclude the feasibility analysis the end-point constraints of the 
problem must also be checked (in addition to the visual check) in order for a 
solution to be feasible. A sample of a feasible solution is provided in Figure 1 
and a sample of an infeasible solution is provided in Figure 2. 





35 |. M. Ross and F. Fahroo, “User’s Manual for DIDO 2001.” 
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Figure 1: | Sample Feasible DIDO Solution 
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Figure 2: | Sample Infeasible DIDO Solution 
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lll. ONE-TO-ONE TRANSFERS 


A. INTRODUCTION 


The purpose of this work is to select the optimal trajectory to rendezvous 
with an asteroid, given a group of asteroids. The selected solution must be 
verified as the optimal. In order to verify the selection, it is necessary to 
determine the cost for the spacecraft’s one-to-one rendezvous with each 
asteroid. These costs can be compared, the minimum determined, and this 
minimum value becomes the basis for optimality for the selection verification. 
This section summarizes the results of the one-to-one spacecraft-asteroid 


rendezvous. 
B. ONE-TO-ONE RENDEZVOUS RESULTS 


The first step taken in this work was to determine the cost for a one-to-one 
rendezvous between the earth and a single asteroid. This analysis is conducted 
for a reasonable number of asteroids (five) for comparison. The selection criteria 
would then be introduced over this group of five asteroids such that the selected 
asteroid and subsequent trajectory could be compared to the baseline 
established with the one-to-one transfers. 


To execute a one-to-one rendezvous, the model of the system simply 
includes the Sun, the Earth, the spacecraft and the particular asteroid in 
question. The dynamics of these bodies is inputted into DIDO to produce a 
solution based on the boundary conditions specified. An initial test was run with 
15 nodes on the final five asteroids in the data set. These five asteroids were 
chosen without any specific intent or special interest. The results of this test are 
shown in Table 4. The cost is displayed in canonical units (mass units per time 
unit). From the tabulated values, it is clear that the minimum cost for one-to-one 


transfer occurs between Earth and asteroid number 910. 
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C. VERIFICATION OF THE RESULTS 


As Feasibility Analysis 


The feasibility plots for these one-to-one rendezvous are shown in Figure 
3 through Figure 7. The 15 node solution initially generated was bootstrapped as 
a guess to create a 60 node solution with greater accuracy. The results of the 60 


node solutions are shown in Table 5. 


Table 4: One-to-One Transfer Costs and Visible Feasibility 




















Asteroid Number | One-to-one Cost | Visibly Feasible 
906 0.15492 Yes 
907 0.20091 Yes 
908 0.094657 Yes 
909 0.14122 Yes 
910 0.08575 Yes 

















The minimum cost is once again found to be in a_ trajectory 
rendezvousing with asteroid 910. These costs have a higher level of confidence 
as they are obtained using more nodes. The increase accuracy is evident in the 


feasibility plots shown in Figure 8 through Figure 12. 
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Table 5: Bootstrapped One-to-one Transfer Costs and Visible 




















Feasibility 
Asteroid Number | One-to-one Cost (MU/TU) | Visibly Feasible 
906 0.1243 Yes 
907 0.26218 Yes 
908 0.16541 Yes 
909 0.14527 Yes 
910 0.067906 Yes 

















In order to verify that the plots of the states in Figure 3 through Figure 12 were 
generated based on the proper propagated solutions, the control variables were 
propagated and checked for their feasibility as well. This confirms that the states 
are being generated by the same control solution that DIDO generated. Although 
the 15 node solutions were found to be visibly feasible, there was room for 
improvement. This is evident from the plots of the feasibility of the control 
solution. In some of the plots, the propagated solution passes in the same 
general shape that the nodes create, but it does not directly pass through each 
one. This shows a lack in the number of nodes used to find the control solution. 
As such, a higher fidelity solution was necessary. The 60 node solution, 
bootstrapped from the 15 node solution, shows a much higher level of 
confidence. In Figure 8 through Figure 12, the propagated solution passes 
directly through every node, with high accuracy and little room for variation. 
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Figure 3: —_‘- Visual Feasibility Plot for Asteroid 906 with 15 Nodes 
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Figure 6: Visual Feasibility Plot for Asteroid 909 with 15 Nodes 
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Figure 9: ‘Visual Feasibility Plot for Asteroid 907 with 60 Nodes 
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Figure 10: Visual Feasibility Plot for Asteroid 908 with 60 Nodes 
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Figure 11: Visual Feasibility Plot for Asteroid 909 with 60 Nodes 
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Figure 12: Visual Feasibility Plot for Asteroid 910 with 60 Nodes 
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Although Figure 8 through Figure 12 imply that the obtained solutions are 
visually feasible, more analysis must be conducted before the solutions can be 
confirmed as feasible. The solution must meet all of the boundary conditions set 
in the problem formulation (i.e. the spacecraft trajectory must start at Earth and 
end at the asteroid). 


2. Checking the Boundary Conditions 


DIDO outputs allow the user to easily check if the boundary constraints 
are met by the solution. The solution variables are all available and can be 
inserted into the events file to check if the end-points constraints are satisfied. 
The value of each constraint must fall between the limits established in the main 
file. There are seven of these end-point constraints: initial position, velocity, and 
mass of the spacecraft, final position and velocity of the spacecraft, mean 
anomaly of the Earth, and the asteroid. In order to confirm the feasibility of the 


solution, all of these constraints must be met. 


The solutions generated in Section B do not meet the constraints on the 
final position and velocity of the spacecraft. These values are shown in Table 6. 
Table 6 displays the boundary conditions as they appear in the events file. The 
initial position takes the difference between the position of the spacecraft and the 
Earth at the first node. 


(Xie, - Xeart ) 5s War a Y ert ) + (Zar = Zea ) (3.1) 
Initial velocity is the difference between the velocities of the spacecraft and the 
Earth at the first node. This value has a range because the spacecraft has an 


initial escape velocity. 


2 2 2 2 2 2 
(vi... ~ Yi! a (vi... ~ VS raias (v2... = Meee (3.2) 


Spacecraft mass is the difference between the initial spacecraft mass and the 
initial mass generated by DIDO. 


M, — Mp 


(3.3) 


calculated 
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Final position is the difference between the position of the spacecraft and the 
asteroid at the final node. This value has a range to allow for a rendezvous 


tolerance. 

ar = Lavy ) a (Ver = Vast ) + (2566 = Zag ) (3.4) 
Final velocity is the difference between the velocity of the spacecraft and the 
asteroid at the final node. This value also has a range to allow for a rendezvous 


2 2 2 2 2 2 
(Vi. a Vp = (Vi. 2 eiah ‘3 Ca Veg 


Mearth refers to the mean anomaly of the Earth. It is calculated as the difference 


tolerance. 


(3.5) 


between the mean anomaly based on the given data and the calculated mean 
anomaly from DIDO all at the first node. 
M 


Earth E, + €earth sin(E, ) (3.6) 
Similarly, Mast refers to the mean anomaly of the asteroid and is calculated in the 


same way, only at the final node. 























Mag — =; + Cage ~ SiN(E; ) (3.7) 
Table 6: Boundary Conditions for the 60 Node Solutions 
Aces Initial Initial Spacecraft Final iby Messr Mac 
Postion Velocity Mass Position Velocity 

Lower Limit 0 0 0 0 0 0 0 

Upper Limit 0 122.5 0 4.47E-13 1.00E-10 0 0 
906 8.84E-7 122.5002 0 0.0012 7.47E-6 -3.99E-14 | -2.84E-9 
907 5.54E-6 | 122.5000 | 1.87E-13 0.0013 3.75E-4 5.51E-13 | -4.33E-7 
908 5.94E-7 | 122.5000 0 2.06E-4 3.24E-6 -1.15E-12 | 1.49E-10 
909 3.63E-6 | 122.5000 0 1.48E-4 8.13E-6 -1.55E-11 | 1.58E-11 
910 5.29E-7 122.5000 | 4.89E-15 8.75E-4 7.32E-4 3.30E-14 | 6.66E-14 
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It is evident from this table that not all of the end-point constraints are met 
by these solutions. Many of the values are so small that the truncation and 
round-off errors inherent in numerical techniques are the cause of the variance. 
However, when the final position and velocity are observed more closely, there is 
an obvious disparity between the limit values and the values generated by the 


solution. 


This disparity can also be demonstrated graphically by plotting the 
trajectory of the spacecraft and the orbits of the Earth and the asteroid. These 
plots are shown in Figure 13 through Figure 17 for the one-to-one rendezvous 
with asteroids 906 through 910 respectively. These figures all confirm that the 
spacecraft trajectory is not ending at the asteroid orbit. Thus, the solutions 
produced in Section B are confirmed to be infeasible. There is no need for an 


optimality analysis of infeasible solutions. 
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Figure 13: 3-D Plot of Spacecraft Trajectory with respect to the Earth and 
Asteroid 906 Orbits 
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Figure 16: 3-D Plot of Spacecraft Trajectory with respect to the Earth and 
Asteroid 909 Orbits 
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Figure 17: 3-D Plot of Spacecraft Trajectory with respect to the Earth and 
Asteroid 910 Orbits 
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D. CORRECTING THE MODEL 


As a result of the analysis presented in Section C of this chapter, the 
solutions for the one-to-one spacecraft-asteroid rendezvous do not pass the 
feasibility and optimality tests. There must, therefore, be a problem with the 
model as it was coded. This section focuses on identifying the error and the 
steps taken to correct this error. 


1. Problem Identification 


When the model is suspected of containing errors, it is necessary to check 
every line of the code for typing and/or mathematical errors. If this line by line 
check fails to unearth the inaccuracies of the model, the dynamics and events 
files should be checked for accuracy. When the events file that was used to 
generate the solution in Section B was reviewed it became evident that the two 
variables which were falling outside of the bounding values (i.e. final position and 
velocity) were not scaled correctly. This is also evident in Table 6 as the 
numbers enforcing these constraints were on the order of 10° and 10°° 
respectively. These values are so small they are not within the DIDO 
computation tolerance. In addition, any round-off or truncation error could cause 
a severe alteration in the outcome. These endpoint conditions must therefore be 
scaled such that the numbers enforcing them can produce a working solution. 


2. Problem Correction 


Having identified a possible source of the error in the model, the 
necessary steps to correct the problem were studied. The problem was re- 
scaled so that all of the constraints are working within reasonably scaled 
numbers, so as to avoid errors caused by round-off or truncation errors. 
Fortunately, only the two constraints pertaining to the two final values in question 
need to be scaled, as the other values are already within reasonable orders of 
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magnitude. These two constraint equations were, therefore, entirely multiplied by 
a common multiplicative factor which would raise their order of magnitude. Both 
the constraints desired bounds and the functions which generate the constraint 
values must be multiplied by the same factor. When properly executed, the 
orders of magnitude of these two constraints are raised, without altering the 
configuration of the problem and without a need to completely identify new 


canonical units. 


This technique for attempting to correct the model’s deficiencies was 
utilized for the same five asteroids. The scale factors were adjusted to obtain a 
solution that meets all of the end-point criteria. Based on the orders of 
magnitude before this constraint scaling, a good starting point could be 
established for bringing the numbers to a reasonable magnitude. Eventually, 
using a scaling factor of 1x10° for the final position value and a scaling factor of 
1x10* for the final velocity value, we obtained one-to-one rendezvous trajectories 
to asteroids 906 through 910 that meet all of the end-point criteria and have a 
feasible solution. The feasibility plots are shown in Figure 18 through Figure 22 
and the resulting costs are given in Table 7. These solutions provide trajectories 
which hold the end-point values within the constraints. These values are shown 
in Table 8. 

















Table 7: Rescaled One-to-One Rendezvous 
Asteroid Number | One-to-one Cost | Visually Feasible 
906 0.1497 Yes 
907 0.2115 Yes 
908 0.1292 Yes 
909 0.1369 Yes 
910 0.0915 Yes 
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Table 8: 


Boundary Conditions for the Rescaled One-to-One 





























Rendezvous 
ee ceeke Initial Initial Spacecraft Final Final eee. Mac 
Postion Velocity Mass Position Velocity 

Lower Limit 0 0 0 0 0 0 0 

Upper Limit 0 122.5 0 0.01 0.01 0 0 
906 6.27E-05 | 122.4954 0 0.0047 0.01005 5.09E-11 | -3.09E-11 
907 1.65E-07 122.5 0 0.0068 0.01003 4.50E-13 | 4.33E-13 
908 4.64E-07 122.5 4.80E-08 0.01002 0.0085 -7.54E-14 | 4.99E-12 
909 2.53E-06 | 122.4994 0 0.0099 0.0098 -1.25E-11 | -1.06E-13 
910 8.46E-06 122.5 -2.20E-13 0.0089 0.01 1.50E-11 2.26E-11 


























The plots of the trajectories show the spacecraft rendezvousing with the 


asteroid as desired. This is shown in Figure 23 through Figure 27. The data 


shown is only for the 30 node solution bootstrapped from an initial solution of 15 


nodes. These solutions meet the conditions for feasibility. 
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Figure 18: Rescaled 30 Node Visual Feasibility Plot for Asteroid 906 
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Figure 19: Rescaled 30 Node Visual Feasibility Plot for Asteroid 907 
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Figure 20: Rescaled 30 Node Visual Feasibility Plot for Asteroid 908 
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Figure 21: Rescaled 30 Node Visual Feasibility Plot for Asteroid 909 
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Figure 22: Rescaled 30 Node Visual Feasibility Plot for Asteroid 910 
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Figure 24: Rescaled 3-D Plot of Spacecraft Trajectory from Earth to 
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Figure 25: Rescaled 3-D Plot of Spacecraft Trajectory from Earth to 


Asteroid 908 





Asteroid Orbit 

Earth Orbit 

Sun 

* Spacecraft Trajectory 














Final Spacecraft Position 





Rescaled 3-D Plot of Spacecraft Trajectory from Earth to 
Asteroid 909 


52 


Figure 26: 





Asteroid Orbit 
Earth Orbit 
*K Sun 
* Spacecraft Trajectory 














Final Spacecraft Position 





Initial Spacecraft Position 





Y (AU) aa ae 


Rescaled 3-D Plot of Spacecraft Trajectory from Earth to 


Figure 27: 
Asteroid 910 


Now that the model is producing feasible solutions, the work is focused on 
producing a selection routine that will select the optimal asteroid with which to 


rendezvous. 
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IV. ASTEROID SELECTION 


A. SELECTION LOGIC 


With a base of knowledge for the expected solution, based on the data 
generated using the one-to-one Earth-asteroid rendezvous, the problem can now 
be formulated to include multiple asteroids. The desired outcome is to find one 
asteroid from the group with the least cost. The selection logic is implemented 
into the events file. A simple selection algorithm can be created that will 
automatically select an asteroid by finding the states of the asteroids at the final 


time, the final positions (x,,y,,Z,) and velocities (v, ,v,,v, ) of the asteroids are 


Xp? We Ze 
known. Comparing these values to the final position and velocity of the 
spacecraft, the asteroid with the matching final state values can be determined 


and selected as the target asteroid. 


. 2 2 2 2 2 2 
ney ( Vets V iss - Ge Yee = (vi. Vou )| (3.8) 


This is a very simple algorithm designed to fulfill the asteroid selection 


requirement of the problem. 
B. SELECTION RESULTS 


Using this selection method, DIDO is now run to find an optimal solution to 
any one of the five asteroids. Knowing the cost to go to each individual asteroid 
(Table 7), it is expected that DIDO will provide a control solution for a rendezvous 
with asteroid 910 based on the results of the one-to-one transfers. The results, 


displayed in Table 9, show that asteroid 910 was not selected. 


Table 9: DIDO Output for a Model with Five Asteroids 





Asteroid Chosen Cost Visually Feasible 





907 0.34364 No 
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This is a very surprising result. The solution was found to be infeasible based on 
the visual check (see Figure 29). The cost of this solution is much greater than 
expected. Also, DIDO selected a different asteroid than expected. This solution 
is obviously not the optimal trajectory for an asteroid rendezvous. This selection 
technique did not produce a working solution. A new approach is, therefore, 


necessary. 


The optimal trajectory is most likely the trajectory which requires the least 
Av for the spacecraft, as it would require the least fuel and is closest to the 
Earth. The velocity of the asteroids at the final time can be compared to the 
velocity of the Earth at the initial time to determine the asteroid requiring the least 
Av for rendezvous. 
min( V5, Veae) (3.9) 
This is another technique for asteroid selection. This method of selection takes 
into account only the mass optimization facet of the problem at hand, however. 
DIDO is run again, implementing this new selection process; the expected 
selection is still Asteroid 910. This technique chose asteroid 907 again. The 
data is shown in Table 10. Unlike with the technique used earlier, however, this 
solution is feasible. The feasibility plot is shown in Figure 30 and the endpoint 
constraints in Table 11 and Figure 28. 


Table 10: DIDO Output for a Model with Five Asteroids with new 
Selection Technique 





Asteroid Chosen Cost Visually Feasible 


907 0.24966 Yes 
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Table 11: Endpoint Constraints for new Selection Technique 














; Initial Initial Spacecraft Final Final 
Asteroid ; ; ae ; Mearth Mast 
Postion Velocity Mass Position Velocity 
Lower 
0 0 0 0 0 0 0 
Limit 
Upper 
0 122.5 0 0.01 0.01 0 0 
Limit 
907 2.299E-7 | 122.5000 | 2.265E-14 | 0.01005 0.01007 | 1.20E-14 | 2.77E-13 
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Figure 28: 3-D Plot of Spacecraft Trajectory with respect to the Earth and 
Asteroid 907 based on the New Selection Technique 


While DIDO produced a feasible solution, the cost of the solution is much 
greater than desired. The one-to-one rendezvous with Asteroid 907 had the 
highest cost of all the one-to-one transfers. Either the selection criterion is 
flawed or DIDO has located a local minimum. In order to further understand this 


result, a discussion on global optimization is necessary. 
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Figure 29: Visual Feasibility of Solution with Selection of Asteroid 907 
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Figure 30: Feasibility of Selected Asteroid Utilizing New Selection 
Technique 
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C. GLOBAL OPTIMALITY 


The purpose of this section is to familiarize the reader with the difference 
between a local and global optimum and how this was utilized in this work. A 
function may have many local minima; however, they can have only one globally 


minimum value. This concept is shown in Figure 31. 





Local Minimum 


Local Minimum 





Local Minimum 


Global Minimum 

















Figure 31: Demonstration of Global Minimum 
It is clear that this illustrative function in Figure 31 has five maxima and four 
minima. Except for the end-points, these points are where the slope of the 


function is zero Y =0), The smallest of the minima is known as the global 
ix 


minimum, it is the smallest value which the function takes over the entire interval. 


It is this global minimum that we are searching for in this work. 
1. DIDO and Global Minima 


DIDO searches for an optimal solution. However, once it has finds a local 
extremal, it will stop its search and output the result. If the function under 
investigation has several local minima, as in the function shown in Figure 371, it is 


possible that DIDO would find only a local minimum. The version of DIDO used 
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in this work does not search for the global optimum. As a result, if the user wants 
to find a global optimum, it becomes necessary to find a technique that would 
lead to this value. This is done by repeatedly running DIDO while introducing 
another constraint: an upper bound on the final cost that is lower than that the 


previously calculated cost. 


Earlier, while discussing the model used to find the solution of the 
problem, the events file was described as containing the specific boundary 
constraints of the problem. To help DIDO find the global optimum, and confirm 
global optimality, the final cost of the maneuver is added as an additional 
constraint and the problem is resolved. If we repeatedly constrain the values of 
the cost function to less that what was found in previous iterations until there is 
no longer any feasible solution the cost of the prior run was the global minimum 
and therefore the global optimal trajectory is determined. Using this technique, 
DIDO has to be run several times, once for each new decrease in the bound on 
the cost. This can become a time consuming practice, but its payoff is the 
acquisition of the globally optimal solution. This technique for searching for the 
global optimum will be called gradual cost-constrained optimization. 


2. Gradual Cost-Constrained Optimization Results 


The results found during the initial testing of the selection methods were 
discouraging. We can use the result of the selection as a jumping off point to 
start the search for the global minimum using gradual cost guiding. We know 
that there a minimum exists which is less than 0.25 based on the data collected 
during the one-to-one transfer study. By forcing the cost function to be less than 
0.1245, DIDO’s output must be a trajectory to asteroid 910 or a different asteroid 
with a lower cost than the one-to-one rendezvous. From this we should get a 
result to continue our search with. Once DIDO no longer finds a feasible result 
within the constraints, we know that we had attained the globally optimal 
trajectory in the previous run. This value should be equal to or less than the 
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minimum cost found using the one-to-one transfer performed earlier. The results 


from the tests using gradual cost-constrained technique are shown in Table 12. 


Table 12: Selection Among Five Asteroids (906-910) with the Gradual 
Cost-Constrained Technique 





Limit on Cost | Asteroid Chosen Cost Visually Feasible 











No Limit 907 0.2499 Yes 
0.125 909 0.0911 Yes 
0.09 907 0.2618 No 




















It can be seen that from the initial starting point, the introduction of the 
upper bound greatly reduced the cost. The cost produced by this trial is slightly 
less than that generated by the one-to-one transfer to asteroid 910, the lowest 
cost previously found. When this cost was used as the upper bound, DIDO 
produced a solution with a cost greater than the bound, immediately nullifying the 
solution by being outside of the end-point constraints. The solution generated on 
the previous run is therefore suspected to be the optimal control because no 
solution could be found with a lower cost. 


D. VALIDATION OF RESULT 


Before concluding that the control solution produced in the previous 
section is a globally optimal trajectory it still must undergo a through verification 
process, the same way we did for the one-to-one rendezvous in Chapter Ill. The 
validity of the solution will be examined from several angles. First, the solution 
must meet all of the boundary conditions set up by the problem. In addition, the 
propagated control should show the feasibility of the results. It also must meet 
the necessary conditions for optimality, as developed in Chapter II. 
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Using the same techniques as before, the validity of this solution must 
also be checked for visual feasibility. The result generated earlier, in Section 
C.2, rendezvousing with asteroid 909 resulting in a cost of 0.08564, is shown to 
be visually feasible based on Figure 32. In addition, the boundary conditions 


must be met. These values, shown in Table 13, demonstrate that the end-point 


conditions are not met. 

















Table 13: | Boundary Conditions for the Selected Trajectory Solution to 
Asteroid 909 
Initial Initial Spacecraft Final Final 
Asteroid ; ; Mearth Mast 
Postion Velocity Mass Position Velocity 
Lower Limit 0 0 0 0 0 0 0 
Upper Limit 0 122.5 0 0.01 0.01 0 0 
909 1.386E-4 | 122.1785 0 3.5556 0.1604 -4.8373E-9 -2.3444 





























The initial velocity, position, mass, and mean anomalies were within the 
bounds of the problem. The final position and velocity, on the other hand, were 
not within the necessary bounding values. This problem is not evident, however, 
in the plot of the trajectory of the flight with respect to the orbits of the earth and 
asteroid as the spacecraft ends on the asteroid’s orbit. This plot is shown in 
Figure 33. 


indeed start at earth, but it does not rendezvous with the asteroid as the final 


It is evident from the values of Table 13 that the spacecraft does 


point on the trajectory does not lie on the orbit of the asteroid. This solution is, 
therefore, infeasible. Checking the solution against the necessary conditions is 
not necessary as the solution has been found to be infeasible. This selection 
technique has been proved to be inadequate as the optimal solution was not 


generated through its use. 
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Figure 32: Visual Feasibility of Selected Trajectory Solution to Asteroid 
909 
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Figure 33: 3-D Plot of Spacecraft Trajectory with respect to the Orbits of 
the Earth and Target Asteroid 909 
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V. CONCLUSIONS 


A. INTRODUCTION 


The purpose of this chapter is to give insight gained from working on this 
project as well as guidance for the continuation of work on this topic. There is a 
great deal yet to be done with this problem as it has not yet been fully solved. 
Hopefully the reader can gain some advantage when undertaking the 
continuation of this work based on what was learned while this research was 


being conducted. 
B. SUGGESTIONS 


The first and most important step for anybody trying to do work with DIDO 
is to familiarize themselves thoroughly with the way that it works. DIDO is a very 
powerful tool for finding optimal control solutions, if it is used properly. DIDO 
requires the problem to be formulated in a specific way in order to get the correct 
output. Having a strong understanding of the requirements and structure used 
by DIDO will relieve the researcher of a great deal of headache and time lost due 
to lack of understanding and errors. For more information regarding the use of 
DIDO, see the DIDO user’s manual (Reference 8). 


If not properly scaled, DIDO can take a very long time to run. Some of the 
trials completed for this work required over ten hours to run. These were run on 
a single PC without any special computing power. The runtime may be 
decreased if a more powerful computing source with several cores were used to 
run the iterations. Another way to decrease runtime is through improved scaling 
and elegant simplifications of the code. Regardless of manner in which it is 
accomplished, this work would be greatly simplified and improved if DIDO’s 
runtime could be decreased. 
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Furthermore, a great deal can be learned from the oversight present in 
this work. The focus, here, was to find a selection process which would choose 
the optimal asteroid with which to rendezvous. The desire to find this value 
overpowered the logical progression of solving the problem. The trajectories 
created in the initial one-to-one Earth-asteroid rendezvous were taken as correct 
simply because the propagation of the control solution provided a visually 
feasible answer. As a result, the error in the model was not found until later. 
Consequentially, a great deal of time was lost once the error was found when 
checking the validity of the selected solution through checking the end-point 
constraints. | Every solution must be checked for feasibility, constraints 
satisfaction, and necessary conditions for optimality before any assumptions of 
their optimality can be made. Any future work should be sure to correct this 
oversight present in this work. 


C. FUTURE WORK 


There is a great deal of work left to do on this project. Once the accuracy 
of the model is confirmed, the next steps can be taken. A selection process must 
be developed that can account for the complex system at hand. The problem 
has multiple layers in that the best asteroid must be selected and the best 
trajectory to that asteroid generated. The techniques used in this research to 
attempt asteroid selection failed to solve this complexity. The first technique 
used assumed linear distances could be used to select an asteroid. The bodies, 
however, move non-linearly as regulated by the equations of motion. The 
second technique, minimizing the Av of the motion, failed to take the time 
optimization component of the problem into account. Other techniques need to 
be explored for the selection process. 


Another source of more progress in this work is the propagation of the 
control. In practice, the only controls available in this problem are T,, Ty, and Tz; 
therefore, T, should be calculated directly from them. When propagating the 


solution, however, the 7 used is the value generated by DIDO. Although the 
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constraint on T? =T, +7/+T? should be enforced when calculating T,, Ty, and 


Tz, it is sometimes not quite as accurate as desired. More accurate knowledge of 
the feasibility of the solution could be gained by propagating the solution with T 


calculated from T,, Ty and Tz instead of the T generated by DIDO. 


It should also be noted that this work focused only on the first step of the 
GTOC2 problem. Once the spacecraft has rendezvoused with one asteroid in a 
group, GTOC2 requires it to rendezvous with three other asteroids in other 
groups. The selection of which group to visit first could be the next step taken 
once this problem has been sufficiently solved. Needless to say, a great deal of 
work still remains to be done in this problem and in the subsequent components 
of GTOC2. Trajectory optimization is an exciting and valuable field, the 
importance of which is growing exponentially. As technology advances, 


mankind’s ability to improve upon the foundations of the past progresses as well. 
D. CONCLUDING REMARKS 


Creating a solution to the Global Trajectory Optimization Competition 
problem proposed by the Jet Propulsion Laboratory proved to be a challenging 
endeavor; one with a great deal of work remaining. It is always a frustrating 
situation when the limitations of time prevent the completion of a project as it has 
here. The foundation for the problem has been established, some tripping points 
have been identified, some obstacles overcome and still much remains to be 


accomplished. 


The goal of this thesis was to develop an optimal trajectory for rendezvous 
with an asteroid in a group. The selection process still requires some attention 
as the optimal asteroid is not being selected based on the tests executed in this 
work. The greatest success of this research was the generation of one-to-one 
Earth-asteroid rendezvous trajectories. Without this data, the asteroid selection 
technique cannot be confirmed. The generation of these trajectories also 


demonstrates that the problem is properly formulated and that DIDO can produce 
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feasible solutions. The gradual cost-constrained optimization technique could be 
a valuable tool for anybody pursuing this problem. The use of this technique is 
not unique to this problem. Any optimal control problem that may encounter 
issues differentiating local and global minima could utilize the gradual cost- 
constrained technique to try to positively identify the global minimum. The 
gradual cost constrained technique can be used throughout the optimal control 
field. Any DIDO user can easily implement this useful tool to enhance their 
solution. Finding a global minimum is, by definition, the goal of every optimal 
control problem. In certain circumstances, this is a very useful method which 


may prove to be useful in the search for this illusive value. 
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ID# 


Group 

2011542 
2001038 
2000624 
2016560 
2001902 
2003552 
2002456 
2015278 
2001173 
3117599 
2001437 
2014569 
2000588 
2002223 
2001754 
2012714 
2007119 
2014268 
2003134 
2000659 
2015436 
2001873 
2012444 
2000225 
3091801 
2001172 
2001208 
2001867 
2006984 
2032511 
2002483 
2002920 
2001345 
2000911 
3079876 
2003709 
2000361 
3046844 
2009799 
2002674 
2002893 


Semi-Major 
Axis 
(AU) 


1 
3.9501468 
3.9619932 
5.2272807 
5.0765482 

3.971589 
4.2298497 
5.1372368 
3.9846457 
5.3227896 
3.0704328 
5.1600747 

3.977245 
5.1947906 
5.1985031 
3.9500648 
5.2062335 
5.2067455 
5.2873272 
3.9793318 
5.1922184 
5.2056216 
5.2491321 
5.2602165 
3.388867 1 
7.2342917 
5.1915116 
5.2374808 
5.133081 1 
3.9692255 
5.0492697 
3.9660655 
5.1130804 
3.9798138 

5.25414 
3.1268516 
5.2609102 
3.9554449 
2.8266851 
5.1942001 
5.1717545 

5.17689 


APPENDIX A: ASTEROID DATA FILE 


Eccentricity 


0.2391642 
0.227483 
0.023528317 
0.0408516 
0.2236554 
0.71337656 
0.074439 
0.2157149 
0.1371044 
0.5119076 
0.0433457 
0.2859047 
0.1465778 
0.0140337 
0.1672029 
0.0356216 
0.1031099 
0.0906027 
0.2201713 
0.1166258 
0.0439938 
0.0922724 
0.0713189 
0.2668163 
0.6863551 
0.10319405 
0.0909268 
0.0437056 
0.1871492 
0.4282629 
0.2763629 
0.026957 1 
0.1808098 
0.065868016 
0.5748146 
0.0614257 
0.21241778 
0.85854077 
0.0479145 
0.0678804 
0.0764874 


Inclination 


(deg) 


6.87574 
9.22988 
18.193628 
15.29915 
12.49318 
30.907935 
13.90266 
9.29657 
6.90993 
21.62245 
20.52178 
10.90696 
10.32134 
15.99057 
12.11822 
9.51657 
19.26732 
14.94484 
7.63747 
4.51942 
16.26675 
21.85938 
30.80052 
20.887897 
24.22652 
16.683194 
33.56575 
26.9088 
16.80098 
8.93717 
4.49881 
21.11915 
11.40178 
21.788152 
20.68275 
19.60369 
12.631655 
20.464243 
30.51326 
1.85443 
14.64719 


LAAN 
(deg) 


16.88982 
58.20488 
342.80533 
100.75566 
59.4253 
350.30192 
327.41517 
344.88781 
283.91796 
355.71658 
315.83789 
346.1705 
316.59634 
220.98833 
163.24455 
298.91469 
285.66811 
284.67115 
257.10507 
350.88872 
253.42915 
197.91284 
213.19945 
197.19179 
144.27642 
247 .39536 
48.55332 
283.69667 
206.29149 
285.89591 
252.1627 
230.96163 
137.5123 
338.02074 
245.06986 
187.17563 
18.960258 
326.68865 
259.54994 
179.86277 
108.76807 
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Argument of 
Periapsis 


(deg) 


48.9603 
307.19412 
185.13006 
157.40487 
267 .89541 
316.99246 

93.63287 
64.13255 
40.01722 
225.09482 
129.75256 

37.5375 

132.48447 

51.37102 
111.93204 
163.16174 
119.08349 
123.62387 
163.00334 
341.45996 
178.28135 
356.47324 

64.17095 

104.2599 
240.32669 
49.477769 
295.69051 
358.70429 
247 .90573 
345.49142 
182.86999 
196.35511 
332.82006 
80.347917 
115.06575 
245.59639 
68.162409 
138.58443 
113.20009 

37.74182 
171.07065 


Mean 
Anomaly 


(deg) 


229.49648 
201.38868 
151.0004 
261.64642 
47.99795 
237.325 
230.61779 
258.62785 
205.92255 
290.06416 
244 .24908 
338.80357 
216.43354 
279.17581 
295.64241 
210.69903 
235.10541 
237.9447 
272.12822 
318.17164 
227 83254 
335.83439 
256.62017 
310.40788 
85.9533 
250.4737 
204.17234 
250.58354 
216.503 
172.48229 
226.77818 
223.31577 
181.8755 
253.5207 
323.76911 
219.05819 
239.42006 
74.835735 
280.79516 
321.3662 
241 .30879 


Epoch 
(JD) 


Asteroid 
Cluster # 


nn nn nn ne ee ee ee ee ee ee 


2011396 
2000617 
2002241 
2003063 
2004834 
3169278 
2002959 
2020898 
2005209 
2002760 
2004489 
2016070 
2002357 
2005264 
2002906 
2003708 
2004833 
2004086 
2003548 
2006090 
3061681 
2001749 
2002759 
2004754 
2005144 
2004791 
3089425 
2004543 
2002363 
2000944 
2005130 
3081550 
2004836 
2004063 
2002797 
2003317 
2011351 
2005012 
2001583 
2005025 
2002260 
2001404 
2001362 
2100004 
2011395 
2002207 
3035962 
2016974 
2004837 
2007641 


5.2054308 
5.2266795 
5.2059683 
5.1854315 
5.2357128 
3.1343661 
3.9440841 
4.2287314 
5.1533221 
3.9834018 
5.2761594 
5.123581 1 
5.1921105 
5.2043321 
3.1616197 

5.220393 
5.2562716 
5.2256622 
5.1604558 
5.3136952 

3.346388 
5.1722902 
5.1700791 
5.1909721 
5.2263273 
5.1792451 
5.3324536 
5.0957247 
5.1655208 
5.7544224 
5.2456965 
3.0440863 
5.1823797 
5.1695887 
5.1057194 
5.2132675 
5.2491115 
5.2652415 
5.1071568 
5.2111349 
5.1896664 
5.3026206 
3.2218943 
2.6090888 
5.2073823 
5.1267457 
4.6900143 
5.2149294 
3.1977504 
5.2186575 


0.0640564 
0.1381751 
0.0670944 
0.0584969 
0.1355501 
0.5204205 
0.2742882 
0.4647309 
0.0495159 
0.1223004 
0.0602848 
0.1236231 
0.0438295 
0.1107648 
0.1137911 
0.1577286 
0.0928909 
0.1209423 
0.0888397 
0.0571189 
0.5398696 
0.1078736 
0.0660091 
0.0080923 
0.2720738 
0.0461581 
0.4300086 
0.0959612 
0.0356465 
0.66019798 
0.0099851 
0.5018147 
0.1077418 
0.1181985 
0.0878189 
0.1259699 
0.0638223 
0.0857824 
0.0522215 
0.0753205 
0.0439744 
0.1135323 
0.3698017 
0.69471344 
0.0672427 
0.0172127 
0.7036003 
0.0702616 
0.1314229 
0.0535495 


12.58719 
22.03587 
16.60687 
12.17397 
28.45078 
20.4243 
5.23244 
45.4942 
9.06609 
13.46141 
22.16378 
16.25994 
2.66967 
33.57519 
30.6904 
13.36738 
34.68256 
21.7387 
8.07572 
20.18055 
22.15021 
6.08844 
21.96815 
12.34371 
8.89879 
25.94351 
42.44729 
14.73139 
32.19648 
42.532955 
15.70813 
16.79417 
19.41902 
18.94705 
22.39234 
2788043 
11.5758 
4.99793 
28.55192 
11.01688 
17.78242 
18.00703 
24.2039 
16.289474 
24.14544 
6.81023 
10.21541 
15.01434 
28.22335 
34.68862 


195.68829 
44.35812 
267.98621 
287 .88023 
76.07509 
252.16872 
121.33837 
293.29012 
322.75689 
352.816 
86.70204 
300.91274 
179.31425 
121.90801 
84.60812 
291.17902 
101.76168 
54.95767 
43.54071 
328.52802 
277.03743 
341.00747 
171.23528 
155.23002 
322.85738 
261.4468 
213.05295 
325.4165 
211.83521 
21.557042 
242.53246 
260.38958 
82.04271 
113.52578 
69.9464 
135.9195 
251.09162 
34.82759 
221.3685 
347.86522 
86.57157 
332.97289 
121.39505 
77.431206 
213.23624 
159.1796 
183.5988 
241.62526 
327.27954 
242.05333 
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175.96375 
307.5839 
291.17519 
204.89888 
350.67143 
250.66259 
284 .68371 
234.87212 
104.50996 
155.94779 
5.3875 
352.61268 
72.48815 
358.9272 
294 .92209 
56.85664 
278.74281 
356.03329 
26.8289 
72.38497 
145.21595 
111.19397 
8.41746 
213.07335 
330.63494 
165.16477 
193.03549 
84.13779 
51.98333 
56.72622 
103.5806 
38.40177 
34.03318 
317.49315 
47.69826 
149.13296 
159.70997 
333.0567 
186.33737 
72.09646 
319.93149 
59.33793 
262.11257 
11.737654 
117.37947 
300.66192 
111.51786 
134.19018 
42.92711 
227 .13421 


289.49289 
183.22681 
324.5333 
178.85982 
223.76284 
169.20267 
297 .55373 
212.42781 
257 .03537 
327.81918 
190.45873 
243.85689 
292.0237 
178.74752 
165.70717 
189.19123 
271.01003 
270.21371 
246.50665 
259.09504 
341.7947 
197.09829 
126.14053 
186.6018 
227 .10228 
99.89186 
136.63782 
262.49894 
290.96291 
43.435582 
169.16343 
356.90507 
195.75845 
233.457 24 
175.56344 
258.45877 
276.73337 
293.9424 
260.04144 
233.10794 
252.91054 
267 .83412 
18.88301 
142.26727 
328.72511 
81.0411 
202.03203 
272.57798 
296 .90452 
184.77791 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


SN ne ee ee ey 


2003793 
2004035 
2005259 
2004709 
2009661 


Group 

2000010 
2000016 
2000021 
2000022 
2000024 
2000031 
2000034 
2000035 
2000036 
2000038 
2000041 
2000047 
2000054 
2000055 
2000058 
2000066 
2000069 
2000070 
2000074 
2000075 
2000078 
2000081 
2000086 
2000090 
2000095 
2000097 
2000099 
2000104 
2000105 
2000107 
2000110 
2000111 
2000120 
2000121 
2000125 
2000128 
2000129 
2000132 
2000134 
2000135 
2000136 
2000137 
2000143 


5.1948389 
5.2795635 
5.1886873 
5.1962516 
3.9462035 


2 
3.1366114 
2.9197742 
2.4350514 
2.9091742 
3.1307722 
3.1498399 
2.6851564 
2.9897564 
2.7470479 
2.7424397 
2.7654349 
2.8776367 
2.7122777 
2.7585286 
2.6998953 

2.6451 
2.980066 
2.6162373 
2.7784206 
2.673238 
2.6201204 
2.8537425 
3.1156178 
3.1571193 
3.067408 
2.668392 
2.6637091 
3.1549271 
2.3739092 
3.4773045 
2.7340757 
2.594318 
3.1157251 
3.4573207 
2.7420783 
2.7489736 
2.8677223 
2.6093823 
2.5632069 
2.4280968 
2.2865589 
3.1184487 
2.7624553 


0.0891835 
0.0562528 
0.0731641 
0.0204462 
0.2336894 


0.11799269 
0.13948318 
0.16374276 
0.10280545 
0.13212508 
0.22591347 
0.1087005 
0.228466 
0.30349954 
0.15148566 
0.27184993 
0.13528824 
0.1964305 
0.14476162 


0.043183319 


0.1733656 
0.16704444 
0.1812037 
0.23987342 
0.3047786 
0.20736367 
0.2108515 
0.20751135 
0.15623753 
0.1489283 
0.25700227 
0.19688701 
0.15268145 
0.17644901 


0.078537037 


0.0781141 
0.10084142 


0.059447364 


0.13841548 
0.08125 
0.12733828 
0.21285798 
0.3881 
0.1166245 
0.20653464 
0.0847066 
0.21988728 
0.0703414 


20.92981 
12.134 
15.93135 
25.51675 
12.9873 


3.8423915 
3.095617 1 
3.0644634 
13.710708 
0.7597936 
26.31623 
5.50356 
7.93736 
18.431224 
6.9542317 
15.764855 
4.984635 
11.80389 
7.1847686 
5.0578201 
3.04729 
8.5813425 
11.58479 
4.0751328 
5.00253 
8.6876426 
7.81207 
4.8208178 
2.2197774 
12.99847 
11.78307 
13.858287 
2.7914565 
21.460972 
10.047151 
5.97379 
4.923668 
6.9548919 
7.599757 
4.65594 
6.2542741 
12.218021 
25.05101 
11.588804 
2.3056068 
9.56953 
13.42254 
11.46947 


200.52401 
233.73474 
67.46291 
253.23497 
56.85707 


283.45768 
150.34418 
80.91394 
66.236811 
36.007074 
31.239635 
184.534 
353.8186 
358.47404 
295.91795 
178.16309 
3.2446812 
313.45001 
10.539693 
161.29511 
7.66789 
185.1207 
47.80504 
197.31394 
359.48157 
333.58458 
1.50533 
86.451459 
70.234113 
243.14905 
159.77806 
41.678878 
41.872068 
188.35744 
173.13537 
56.99462 
305.89396 
341.51217 
73.217451 
169.16097 
76.459095 
136.44012 
258.91641 
346.2139 
343.89794 
186.53627 
202.45635 
333.23845 


73 


262.18506 
197.12702 
199.37861 

89.40599 
288.56144 


313.00295 
227.86709 
250.10776 
356.08934 
107.94133 
62.002727 
330.0982 
213.97512 
47.141652 
168.63723 
46.221564 
314.60266 
345.58093 
4.2608489 
34.44369 
43.73999 
289.98554 
255.87084 
174.52383 
339.5614 
151.41908 
50.1751 
307.76761 
242.52717 
155.04622 
268.6762 
196.02815 
30.942654 
56.502696 
309.88106 
281.76827 
165.92105 
231.93546 
295.98151 
110.22377 
302.41385 
108.17944 
254.37088 
83.692584 
340.03522 
132.61377 
106.78731 
250.94068 


206.49881 
232.16471 
37.88812 
213.85689 
8.57451 


56.250366 
101.78424 
231.18992 
3.206616 
257.16605 
356.86228 
178.03603 
58.38133 
25.346092 
86.025664 
226.08638 
204.79556 
81.759176 
93.446378 
352.92475 
46.22178 
137.4484 
78.00096 
15.554851 
3.77453 
330.57112 
129.18942 
247 .06828 
58.626905 
308.59388 
351.69754 
281.54555 
188.30844 
7.3945065 
346.53934 
284.77871 
61.815975 
132.32475 
233.08272 
246.10837 
159.53128 
90.343092 
276.89429 
86.432954 
220.0661 
346.10322 
297 .02006 
174.63933 


54000 
54000 
54000 
54000 
54000 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


A a 


NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNDND LYM 


2000144 
2000145 
2000146 
2000147 
2000156 
2000159 
2000161 
2000163 
2000168 
2000171 
2000173 
2000175 
2000185 
2000187 
2000194 
2000195 
2000200 
2000201 
2000205 
2000206 
2000207 
2000209 
2000211 
2000216 
2000224 
2000232 
2000238 
2000240 
2000250 
2000266 
2000304 
2000313 
2000325 
2000326 
2000329 
2000334 
2000338 
2000342 
2000344 
2000345 
2000347 
2000350 
2000356 
2000360 
2000369 
2000373 
2000375 
2000380 
2000381 
2000382 


2.6544514 
2.6731693 
2.7189094 
3.1342865 
2.7331251 
3.1010391 
2.3795391 

2.367393 
3.3757402 
3.1334027 
2.7416607 
3.1847599 
2.7394441 
2.7321898 
2.6181429 
2.8803618 
2.7369494 
2.6787336 
2.7767208 
2.7403923 
2.2836787 
3.1445726 
3.0409546 

2.797122 
2.6444588 
2.5502931 
2.9086285 
2.6643756 

3.153134 
2.8043578 
2.4030494 
2.3756945 

3.202845 
2.3176268 
2.4764225 

3.885343 
2.9120186 
2.568 1538 
2.5948773 
2.3253732 
2.6131906 
3.1112374 
2.7559652 
2.9986233 
2.6491767 
3.1131734 
3.1231061 
2.6788222 
3.2204713 
3.1156156 


0.23556544 
0.14450786 


0.064905901 


0.0340697 
0.22220581 
0.11097779 
0.1374982 
0.1903735 
0.0673035 
0.12927424 
0.20834644 
0.2324518 
0.1271362 
0.23670797 
0.23614207 
0.04057 
0.13381634 
0.1792631 

0.0361837 


0.040846283 


0.0286959 


0.062655847 


0.16227833 
0.25036608 
0.0460098 
0.1778829 
0.0879187 
0.2068432 
0.1270277 
0.15736045 
0.2214414 
0.1794471 
0.1680242 
0.19054712 
0.0238695 


0.024215343 


0.0201877 
0.128678 
0.31596263 


0.061343146 


0.1641544 
0.1561224 
0.2397877 
0.1817627 
0.0976895 
0.1465925 

0.107686 
0.1132572 
0.0955587 
0.1770362 


4.8084344 
12.636769 
13.074324 
1.93499 
9.7483653 
6.1274027 
9.05318 
4.8058965 
4.63413 
2.5460896 
14.207656 
3.21884 
23.22034 
10.597906 
18.485624 
6.96866 
6.9014922 
5.75756 
10.69477 
3.7803055 
3.80295 
7.1714093 
3.8821037 
13.133755 
5.83838 
6.0708 
12.4025 
2.10451 
12.82149 
13.390532 
15.83975 
11.6459 
8.54318 
23.724144 
15.88431 
4.642853 
6.03807 
7.34603 
18.356075 
9.751305 
11.69434 
24.90012 
8.23048 
11.71228 
12.7077 
15.44736 
15.93372 
6.15571 
12.52422 
7.40193 


76.486474 
77.453037 
84.178218 
248.73039 
242.18207 
134.33009 
18.78424 
160.36447 
206.46806 
100.54761 
148.35779 
21.39951 
153.94657 
21.904126 
159.51944 
7.20269 
324.69997 
157.11692 
212.0672 
145.27928 
29.28743 
0.78538285 
263.75641 
215.66472 
353.0192 
152.52032 
184.19336 
115.22263 
24.05332 
236.00128 
159.22264 
176.82084 
345.28451 
32.334512 
178.54468 
130.2253 
287.63451 
232.74759 
48.25056 
212.79371 
85.84551 
90.19672 
354.86413 
132.65874 
94.38445 
4.0578 
336.67192 
95.21741 
125.34741 
313.60394 
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293.65938 
44.944313 
143.44602 
106.84911 
337.90684 
335.49307 
294 .39304 
297 49222 
168.09446 
58.061999 
228.00507 
321.25122 

224.0974 
195.14152 
162.83374 

123.4987 
86.014728 
181.24502 
171.02616 
302.04137 
192.12819 
250.07104 
174.76297 

179.3494 
284 .05593 

51.12076 
207.19158 
300.50223 

75.74685 
151.25527 

172.3242 
315.81027 

67.86966 
238.50135 

52.19986 
151.07984 
123.54059 

225.4028 
237 38162 
229.68958 

84.86164 

338.9865 

78.80956 
288.55851 
269.55336 
347.77644 
344.59878 
240.54282 
137.52796 
270.78992 


40.157765 
140.69001 
130.23606 
167.74868 
196.55934 
282.20833 
102.63249 
26.693449 
325.53015 
89.493077 
117.48297 
130.21293 
264.43276 
153.41892 
284.62328 
102.20799 
305.98131 
349.42127 
158.14767 
284.70264 

87.89963 

34.32137 
309.57543 
181.60763 
257 .00632 

64.32308 
178.19462 

324.5337 
142.31316 
16.632347 
191.50577 

13.19227 
111.00173 
100.32007 

96.05795 

5.201958 
295.61868 
193.24988 
172.27458 
8.5651683 
248.97015 

249.8531 
333.70592 

27.11377 

10.20211 
202.41227 
192.61115 

21.26849 
102.99105 

81.18543 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNDN NH 


2000386 
2000388 
2000393 
2000395 
2000404 
2000405 
2000407 
2000410 
2000413 
2000414 
2000418 
2000423 
2000429 
2000441 
2000442 
2000444 
2000445 
2000448 
2000449 
2000466 
2000481 
2000488 
2000489 
2000490 
2000494 
2000497 
2000498 
2000508 
2000511 
2000516 
2000521 
2000535 
2000558 
2000559 
2000566 
2000569 
2000583 
2000585 
2000593 
2000602 
2000618 
2000635 
2000654 
2000701 
2000702 
2000712 
2000713 
2000735 
2000751 
2000755 


2.8945485 
3.0067013 
2.7790703 
2.7858859 
2.5943355 
2.5837501 
2.6248052 
2.7285227 

2.584411 
3.5083344 
2.5920701 
3.067624 1 
2.6070974 
2.8060136 
2.3449721 

2.771173 
3.1997043 
3.1381577 
2.5526781 
3.3583767 
2.7389796 
3.1579038 
3.1535619 
3.1685051 
2.9888645 
2.8579034 
2.6503781 
3.1607974 
3.1659016 
2.6796943 
2.7417331 
2.5690206 
2.9064565 
2.7122622 
3.3823202 
2.6561647 
3.1705289 
2.4303409 
2.6975729 
3.0914405 
3.1900659 
3.1443395 
2.2968024 
3.0157244 
3.1943626 
2.5738464 

3.391033 
2.7292584 
2.5502414 
3.1724209 


0.17292904 
0.05887506 
0.33144266 
0.0838058 
0.1983143 
0.24463276 
0.0707861 
0.2365507 
0.3431661 
0.0689192 
0.119635 
0.041046488 
0.12333371 
0.0827663 
0.0714341 
0.17286275 
0.1909852 
0.1843979 
0.1708912 
0.082421992 
0.1582352 
0.1692412 
0.038868229 
0.0989462 
0.0568189 
0.2973859 
0.2251495 
0.014003443 
0.18564446 
0.2735302 
0.2813483 
0.0229292 
0.0431148 
0.064323766 
0.1108082 
0.1828515 
0.1619588 
0.1296983 
0.2179818 
0.2440861 
0.078422036 
0.0786644 
0.23192318 
0.0315655 
0.022868874 
0.1881205 
0.16608003 
0.3214867 
0.153174 
0.1466565 


20.254208 
6.4580588 
14.870607 
3.35143 
14.115702 
11.950979 
7.53482 
10.92226 
18.71569 
9.54215 
6.82571 
11.240443 
9.5273188 
8.1419 
6.06243 
10.278705 
21.37172 
12.71486 
3.08986 
19.16314 
9.8583 
11.500393 
12.976896 
9.26516 
7.07552 
4.82186 
9.50401 
13.36355 
15.938414 
12.95664 
10.59138 
6.78274 
8.3662 
9.3105401 
4.89863 
1.29584 
8.25067 
7.5572 
16.89206 
15.07537 
17.012019 
11.03908 
18.124633 
7.1136 
20.60234 
12.781487 
10.359739 
16.87956 
15.61479 
3.23894 


166.94275 
354.62692 
212.51713 
259.62654 
92.674305 
255.29974 
294.83119 
97.21109 
103.92107 
110.74346 
249.0928 
69.554082 
220.04862 
253.86288 
135.03435 
195.8366 
292.4142 
37.37604 
86.03696 
291.19859 
67.0244 
84.980185 
167.18408 
178.5047 
38.39191 
6.77463 
97.49661 
44.510226 
107.6717 
328.88736 
89.69902 
84.87326 
143.87801 
112.21021 
80.26508 
301.97392 
257 .98239 
180.37753 
76.18743 
331.64478 
111.19981 
183.46648 
278.57048 
244.14352 
289.97403 
231.04829 
217.79282 
43.03 
78.93492 
177.2636 


15 


220.14756 
333.05707 
91.047917 
12.95644 
121.11608 
309.27005 
82.01003 
172.00263 
252.68584 
326.64868 
125.96347 
206.55983 
168.7474 
201.35201 
85.04871 
155.21028 
79.8239 
295.24661 
46.37252 
245 .93322 
348.66307 
68.917033 
12.619788 
196.80207 
216.78947 
2.80643 
241.2688 
179.85303 
338.52068 
258.53574 
316.09633 
68.24631 
317.36518 
130.74 
294.61109 
141.80532 
253.69818 
327.97351 
30.92169 
45.83763 
227 .86543 
219.42795 
214.07158 
313.54037 
352.37054 
181.1623 
136.08384 
310.136 
302.31969 
42.40449 


65.82279 
49.966331 
126.41939 
167.79539 
271.32457 
353.02963 

40.65893 

59.08877 

308.8369 
352.52959 

103.5344 
290.60245 
325.59031 
313.76537 
278.05265 
322.30187 
291 .93822 
107.79207 

342.694 
345.95823 

80.4701 
289.65081 
7.9310261 
120.23447 
177.45085 
228.05837 

75.92796 
129.14648 
247.71727 
296.93286 
246 .32803 

20.05767 

215.4325 
212.48385 
331.04968 

34.01162 
179.04336 

205.455 
286.61962 

302.7661 
239.39252 
168.40125 
112.35919 
228 .69453 
337.72603 
80.890211 
286.71598 
332.75553 
163.86706 
145.37775 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNDND NH 


2000764 
2000772 
2000776 
2000785 
2000786 
2000791 
2000798 
2000814 
2000821 
2000849 
2000860 
2000872 
2000907 
2000931 
2000977 
2001015 
2001028 
2001061 
2001082 
2001093 
2001102 
2001277 
2001445 
2001461 
2001580 
2001606 
2001625 
2001639 
2001794 
2001931 
2001963 
2002379 
2002407 


Group 

2000003 
2000005 
2000006 
2000007 
2000008 
2000009 
2000011 
2000012 
2000014 
2000015 
2000017 
2000018 
2000020 
2000023 
2000025 


3.1850605 
3.0035141 
2.9350442 
2.5692316 
3.1687258 

3.115249 
3.0142313 
3.1498691 
2.7776856 
3.1547846 
2.7945762 
2.7313074 
2.7977888 

3.185335 
3.1153896 

3.208478 
3.3947356 
3.1383114 
3.1216608 
3.1304429 
3.0687473 
2.6992491 
3.1234889 
3.1272276 
2.1968759 
2.6913513 
3.2016194 
2.5719099 
3.1288813 
2.5404735 
2.4215853 
3.1646802 
2.9224974 


3 
2.6676188 
2.5736689 
2.4251131 
2.3854903 
2.2015458 
2.3871957 
2.4522002 

2.334669 
2.5849886 

2.643284 
2.4700724 
2.2955326 
2.4091511 
2.6273733 
2.3999647 


0.1070606 


0.092042448 


0.16086041 
0.2105935 
0.1671064 
0.1996704 
0.0413488 
0.3093096 

0.207103 
0.1953547 
0.1091712 
0.0784015 
0.1637967 
0.2233883 
0.0293926 


0.081469182 


0.1185204 
0.2082506 
0.1808789 
0.2707514 
0.1173315 
0.239092 
0.1772867 
0.0423866 
0.48791122 
0.3150995 
0.2218822 
0.1512801 
0.1552413 
0.27278 
0.21073058 
0.2762638 
0.2211028 


0.25819419 
0.19269048 
0.20172572 
0.23142212 
0.15620731 

0.1214414 
0.10011807 
0.22056537 
0.16811947 
0.18718147 

0.1344116 
0.21870798 
0.14287954 

0.2329636 
0.25544399 


10.07278 
28.783444 
18.245702 

12.73156 

14.55243 

16.38494 

9.22952 
21.83526 
5.37832 
19.48619 
13.31523 
7.36689 

19.57368 

11.44311 

15.19978 
9.4580072 

9.39108 
2.49726 
1.85016 
25.20975 
15.81039 
6.9664 
2.28497 

15.32426 

52.097105 

7.70009 
15.55532 
8.42685 
14.50739 
8.2461 
25.051232 
0.46741 
2.47635 


12.971682 
5.3685634 
14.752027 
5.527283 
5.8884175 
5.5765165 
4.6247051 
8.3638054 
9.1072838 
11.738272 
5.58717 
10.125237 
0.70691821 
10.145257 
21.584123 


259.37313 
64.045057 
79.867084 
72.18331 
89.92703 
130.11717 
214.52891 
88.86011 
209.91865 
228.51247 
309.55072 
194.94024 
43.17347 
111.48085 
75.90784 
120.49107 
63.48771 
91.29504 
148.02332 
55.69696 
216.84653 
247 .25447 
89.30058 
104.84379 
62.32879 
190.77887 
322.144 
324.37971 
221.50437 
182.53508 
106.98322 
151.27404 
342.34897 


170.12215 
141.6853 
138.74385 
259.72283 
110.96399 
68.973442 
125.62803 
235.53818 
86.461325 
293.27326 
125.60807 
150.53455 
206.50811 
67.227831 
214.268 


76 


154.7591 
142.09532 
306.61868 
129.96803 
133.60068 
202.04482 

41.83994 
297.08583 

32.90585 

63.9943 

19.97543 

20.02667 

88.29047 
313.89267 

83.30353 
282.54502 

25.06466 
303.85151 
188.61985 
251.94642 
114.79816 

47.368 
270.99423 
334.56144 
159.50668 
142.30677 
282.75245 
105.46543 
335.55787 
163.92287 
355.82804 
177.62201 

10.80364 


247 .82331 
357.50991 
239.55732 
145.41092 
285.39792 
5.6901086 
195.29555 

69.58648 

96.32871 
97.914615 
136.00231 
228.00031 
255.50551 
59.312781 
90.161994 


100.1314 
205.57057 
249.80628 
128.71164 
115.60687 
217.95885 

8.942 

85.74032 
163.61022 
160.73835 
210.65456 

299.6105 
211.84409 

191.3093 
169.14985 
235.74511 

47.074 
225.10094 

134.3123 
183.03763 
256.27357 
145.46956 
282.75712 
191.21196 
121.45456 
251.45853 
181.27996 

72.49232 
231.28079 

45.11593 
179.40627 
154.27748 

43.83693 


75.986335 
266.07961 
326.18977 
349.83031 
246.60931 
127.8487 
178.97445 
301.90748 
156.36084 
309.03614 
216.015 
43.626263 
346.24964 
351.82574 
33.436221 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


NNMNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNDN NM 


WBWWWWWWWWWWW WW W 


2000026 
2000027 
2000028 
2000029 
2000030 
2000032 
2000033 
2000037 
2000039 
2000040 
2000042 
2000043 
2000057 
2000060 
2000061 
2000063 
2000067 
2000068 
2000071 
2000073 
2000079 
2000080 
2000082 
2000089 
2000100 
2000101 
2000103 
2000113 
2000115 
2000116 
2000118 
2000119 
2000123 
2000124 
2000126 
2000138 
2000149 
2000151 
2000158 
2000167 
2000169 
2000170 
2000172 
2000174 
2000178 
2000179 
2000180 
2000181 
2000182 
2000183 


2.656335 0.086901823 


2.3476948 
2.7780323 
2.5540421 
2.3665367 
2.5879007 

2.865151 
2.6414184 
2.7686996 
2.2678373 
2.4419555 
2.2032648 
3.1493067 

2.393591 
2.9819926 

2.395578 
2.4213289 

2.782613 
2.7549414 
2.6660162 
2.4444158 

2.296534 
2.7598915 
2.5500653 
3.0933126 
2.5829817 
2.7026333 
2.3755464 
2.3806207 
2.7695069 
2.4371741 
2.5812942 
2.6942954 
2.6301198 
2.4389509 
2.4486111 
2.1744038 
2.5917414 
2.8685446 
2.8527952 
2.3578594 
2.5532105 
2.3803732 

2.859052 
2.4599218 

2.970494 
2.7200002 
3.1402839 
2.4179605 
2.7931846 


0.17191778 
0.14824802 


0.072585419 


0.12638171 
0.08299776 

0.3376296 
0.17665525 
0.11419824 


0.046566607 


0.22279671 
0.16794702 
0.1182737 
0.18203332 
0.1678048 
0.12605849 
0.18484185 
0.18536679 
0.176456 
0.041362 
0.19213223 
0.2003401 
0.2244582 
0.18377079 
0.1650934 
0.14151875 
0.0795678 


0.087703611 


0.1915 
0.13755329 
0.1633635 
0.0810231 
0.1218106 


0.076540038 


0.105986 
0.1624635 
0.0653143 
0.0331371 
0.0565242 
0.0336352 

0.13075 
0.0646822 
0.1142941 
0.1458847 
0.0438536 
0.1155056 
0.1690067 
0.1982904 

0.18506679 
0.349925 


3.5621827 
1.5837422 
9.4013833 
6.096443 
2.0974832 
5.5306247 
1.87042 
3.0732741 
10.382913 
4.2556875 
8.5295276 
3.4679469 
15.20008 
3.6020835 
18.21925 
5.7855725 
6.0269995 
7.9716224 
23.255563 
2.37313 
4.6226653 
8.66477 
2.8333141 
16.140829 
6.43021 
10.198982 
5.421 
5.0372639 
11.59731 
3.5690216 
7.74344 
5.77832 
6.42781 
2.9507408 
2.9245 
3.20806 
0.93695 
6.4444 
1.00337 
2.21049 
5.50226 
14.4023 
10.03136 
12.12754 
1.89948 
7.81766 
0.87057 
18.7985 
2.0028901 
26.37466 


45.885083 
94.806111 
144.50295 
356.49859 
307.77461 
220.57549 
8.59314 
7.4126651 
157.17103 
94.292217 
84.397918 
264.93471 
199.33961 
191.80383 
333.7717 
337.91498 
202.72444 
44.182403 
316.10551 
7.23877 
206.80097 
218.82041 
25.636344 
311.64793 
127.33508 
343.4749 
136.27878 
123.59522 
308.99545 
64.03691 
47.7451 
203.73791 
307.95504 
188.18585 
23.47891 
54.95009 
159.64764 
39.04645 
278.55081 
166.44873 
354.82138 
301.46438 
332.0841 
327.80693 
51.19998 
252.12535 
312.70655 
143.59272 
107.27931 
142.01678 


of 


193.16118 
356.77913 
342.55305 
63.459167 
86.72257 
339.79731 
338.24659 
62.695182 
209.57309 
268.90779 
236.6316 
15.950712 
212.88978 
270.41631 
13.81601 
295.63623 
106.29921 
305.38918 
267.4553 
54.66187 
200.35481 
139.11469 
110.38231 
45.00997 
185.87533 
347 .82922 
190.13768 
79.053555 
96.7467 
93.101434 
33.63437 
171.29982 
124.95905 
63.155748 
327.99036 
260.0207 
251.12748 
134.55059 
143.9388 
125.83834 
334.69519 
157.69423 
359.24594 
289.76608 
211.55291 
105.37094 
175.49316 
317.42132 
309.90481 
264.16968 


138.34082 
153.80173 
15.27945 
253.78959 
223.47454 
22.315567 
128.61147 
6.3427201 
36.855297 
220.356 
96.043894 
191.974 
50.32014 
64.514348 
110.25816 
356.80652 
156.01781 
3.378668 
352.31621 
191.9127 
123.7389 
234 .60882 
269.9886 
104.93539 
130.37461 
101.84486 
52.67649 
337.73369 
221 .98288 
318.91078 
198.28777 
338.06631 
16.10407 
230.10546 
91.12117 
143.2047 
253.78577 
7.94691 
3.12319 
28.77731 
350.02999 
300.03171 
288.80348 
316.3647 
281.47495 
104.12798 
231.86403 
113.3454 
156.27731 
253.8036 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


WDBWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW WwW Ww 


2000186 
2000188 
2000189 
2000192 
2000196 
2000197 
2000198 
2000202 
2000204 
2000208 
2000215 
2000218 
2000219 
2000221 
2000228 
2000230 
2000234 
2000235 
2000236 
2000237 
2000243 
2000245 
2000254 
2000258 
2000262 
2000264 
2000270 
2000277 
2000287 
2000288 
2000295 
2000296 
2000305 
2000306 
2000311 
2000312 
2000321 
2000328 
2000339 
2000340 
2000341 
2000346 
2000351 
2000352 
2000354 
2000364 
2000374 
2000376 
2000378 
2000384 


2.3617979 
2.7629704 
2.4499896 
2.404096 
3.1149136 
2.7408956 
2.4593202 
3.0764769 
2.67313 
2.891708 
2.7669402 
2.6657744 
2.3539776 
3.0115033 
2.2013308 
2.3832188 
2.3858256 
2.8829304 
2.8031571 
2.7622571 
2.8611801 
3.1015974 
2.195187 
2.6152135 
2.5523258 
2.802004 
2.1980431 
2.8871937 
2.3527778 
2.7554606 
2.7959558 
2.2287841 
3.1048717 
2.3578262 
2.8978498 
2.7808834 
2.8862686 
3.1079701 
3.0128991 
2.7459333 
2.1992684 
2.7950112 
2.7639952 
2.1939318 
2.8004843 
2.2209348 
2.7796261 
2.2888879 
2.7767367 
2.6511976 


0.1499119 
0.1775287 
0.0370351 
0.24621664 
0.02269419 
0.1602645 
0.2278889 
0.0962809 
0.1719995 
0.0154312 
0.0344195 
0.11716641 
0.2230137 
0.1032757 
0.241223 
0.060890097 
0.2442783 
0.0598554 
0.1873378 
0.0733517 
0.045985314 
0.19769977 
0.1215113 
0.2051338 
0.2141044 
0.1338519 
0.15082857 
0.0873359 
0.023729805 
0.2101927 
0.1703209 
0.1597919 
0.1859286 
0.15049 
0.0078441 
0.1616625 
0.0430939 
0.1126465 
0.09462 
0.1173656 
0.1937417 
0.1025065 
0.156579 
0.1501057 
0.11333869 
0.1490966 
0.0798103 
0.1714106 
0.1295631 
0.1483651 


13.17267 
11.73414 
5.17924 
6.8170145 
7.2609296 
8.79318 
9.30913 
8.82984 
8.27148 
1.74941 
1.69006 
15.226448 
10.84215 
10.8869 
2.53831 
9.4382745 
15.35266 
9.02719 
7.68294 
9.75594 
1.1375042 
5.1773835 
4.51431 
14.29304 
7.70896 
10.43526 
2.3654059 
1.16195 
10.023053 
4.32933 
2.70624 
1.74678 
4.44565 
7.26805 
3.22488 
9.03534 
2.59384 
16.07973 
9.92931 
4.67867 
5.66779 
8.76066 
9.19393 
3.38211 
18.379466 
6.00474 
8.98647 
5.43035 
7.01027 
5.60405 


14.87064 
241.2132 
203.61856 
343.41423 
72.55494 
81.68189 
268.52803 
137.05498 
205.21481 
4.5539 
25.05574 
170.87994 
200.9526 
141.94429 
313.42898 
239.96029 
144.63948 
66.23928 
186.13681 
84.44051 
324.18546 
61.525116 
28.54322 
207.70149 
38.71081 
49.77904 
254.5632 
231.61197 
142.48241 
120.56933 
276.12026 
121.59151 
207.85299 
142.03554 
81.1549 
6.74435 
40.46183 
352.61251 
173.79528 
27.11612 
29.1945 
92.16618 
99.44622 
24742532 
140.45317 
105.61761 
219.23568 
302.25603 
232.7592 
48.10942 
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315.24035 
68.31891 
166.82805 
29.840171 
199.92548 
246 .00807 
88.55903 
0.94421 
55.82695 
125.39114 
321.44768 
60.735184 
142.27416 
195.91626 
18.77463 
139.39921 
192.16925 
211.97564 
174.03894 
201.72651 
108.3985 
327.52007 
233.09692 
154.94955 
24.60432 
340.02132 
80.326577 
135.57413 
120.55732 
83.12444 
148.56188 
252.57079 
260.20059 
167.61862 
40.37368 
260.2729 
30.91652 
101.00197 
160.47059 
42.5169 
293.41228 
289.94435 
31.66235 
144.23938 
7.156191 
312.93786 
27.75383 
316.30919 
156.09891 
33.56521 


18.13794 
295.83763 
133.14401 
221.89748 
346.38828 
318.73523 
290.62198 

72.40597 
137.45271 
191.44426 

34.44622 

77.80577 
336.21897 
104.86267 
353.21101 
230.91656 
134.46103 

170.0847 
258.17166 
197.94377 

28.35885 
193.52253 
329.78941 

38.8324 
119.04846 
217.85622 
210.76351 
264.50517 
223.14593 

123.7304 
226 .83366 
286.30861 
112.25397 
189.61238 
271.97283 

0.17345 
123.78883 

41.63909 
317.36304 
323.31427 

8.27634 
152.33019 
228 .34098 

62.80526 
84.076712 
210.62896 
246 .94554 
318.79015 
167.15252 

73.08764 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
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2000385 
2000387 
2000389 
2000394 
2000397 
2000402 
2000403 
2000416 
2000421 
2000432 
2000443 
2000453 
2000458 
2000459 
2000462 
2000470 
2000471 
2000472 
2000477 
2000478 
2000480 
2000482 
2000487 
2000496 
2000502 
2000509 
2000513 
2000519 
2000529 
2000532 
2000533 
2000534 
2000540 
2000542 
2000548 
2000549 
2000550 
2000556 
2000562 
2000563 
2000565 
2000571 
2000574 
2000579 
2000582 
2000584 
2000599 
2000611 
2000616 
2000619 


2.8472955 
2.7391293 
2.6089236 
2.7601348 
2.6347117 
2.5584357 
2.8107151 
2.7913658 
2.5407221 
2.3691579 
2.2155777 
2.1829344 
2.9945128 
2.6205643 
2.8740387 
2.404738 
2.8861595 
2.5439924 
2.4151108 
3.0152043 
2.6444378 
2.9986012 
2.6712337 
2.1987751 
2.3826518 
3.0645243 
3.0200353 
2.78994 
3.0167164 
2.7705832 
2.9792699 
2.883731 
2.2187693 
2.9068967 
2.2830652 
2.6819981 
2.5885648 
2.465813 
3.0200435 
2.7115251 
2.4441963 
2.4096288 
2.2520296 
3.0098865 
2.6093875 
2.3743584 
2.7703969 
2.9811262 
2.5543067 
2.5202746 


0.12656961 
0.23705267 
0.065160176 
0.229052 
0.2465269 
0.1127408 
0.0966235 
0.21853534 
0.2827883 
0.1462883 
0.040032951 
0.1089484 
0.2423835 
0.2094949 
0.0830898 
0.0933121 
0.23353842 
0.0938108 
0.1882294 
0.0882385 
0.0466305 
0.1030376 
0.086408 1 
0.079568 
0.1791063 
0.0899784 
0.0784204 
0.185968 
0.0952662 
0.17861371 
0.0425945 
0.0570623 
0.0900753 
0.140981 
0.1843779 
0.2607464 
0.2210401 
0.1016516 
0.0947892 
0.2362785 
0.128388 1 
0.2426047 
0.2397801 
0.0827351 
0.2250105 
0.2329791 
0.293821 1 
0.1184348 
0.0578729 
0.0752591 


13.565001 
18.134397 
8.1342143 
6.22412 
12.835543 
11.82138 
9.15501 
12.862244 
7.77214 
12.13154 
4.2304545 
5.5578 
12.62356 
10.29674 
3.1911 
7.22745 
14.98526 
15.79966 
5.28847 
13.17546 
21.29332 
14.46735 
10.23447 
3.78914 
25.00994 
15.41164 
9.71732 
11.01574 
11.02327 
16.313485 
6.55229 
3.27661 
5.57618 
12.06846 
3.87117 
3.96628 
10.11404 
5.23196 
11.12606 
10.24831 
10.99229 
5.22678 
5.68496 
11.02146 
30.01247 
10.7257 
16.67178 
13.44867 
14.96176 
13.78251 


345.24165 
128.31422 
282.55924 
67.37124 
228.26765 
129.5359 
244.84291 
58.208003 
187.53759 
88.87357 
175.55548 
11.82328 
134.95413 
29.57945 
105.32977 
173.30542 
84.095024 
127.26845 
10.73814 
234.02285 
237 .39549 
179.55439 
114.89867 
207.7666 
133.09316 
217.81003 
184.69574 
44.81205 
65.28007 
107.60159 
180.57574 
94.2555 
202.26254 
153.25681 
108.5126 
291.64694 
270.82874 
286.23223 
70.78883 
85.46154 
226.06462 
3.25658 
336.85928 
82.83799 
155.81846 
282.30101 
44.68837 
189.87939 
356.2994 1 
187.56918 
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188.09878 
157.68246 
263.55197 
269.66619 
139.38269 
18.17371 
251.86352 
198.84647 
209.22374 
174.15765 
348.80375 
220.16628 
274.82091 
19.25982 
250.36867 
46.2614 
314.48203 
296.35407 
322.28708 
241 .42639 
212.40356 
88.07899 
281.04144 
258.15505 
19.42627 
157.19742 
222.61276 
303.08502 
336.21389 
76.778716 
40.39599 
333.48245 
337.04541 
214.26388 
320.26734 
156.97231 
44.6313 
177.69124 
261.54824 
336.61987 
290.83268 
27 .65282 
76.87759 
231.62261 
309.97325 
84.784938 
292.98145 
253.66966 
108.68871 
178.39417 


148.57539 
180.64512 
227.94975 
179.60716 
93.681871 
115.8071 
320.86082 
225.9835 
44.4314 
110.98779 
213.66194 
299.0741 
141.13805 
334.1296 
293.74567 
88.91894 
42.91019 
218.67478 
359.09369 
277 .32205 
342.11821 
322.25031 
229.32471 
259.49828 
26.13608 
309.23183 
193.18467 
61.36499 
338.78737 
81.370926 
315.61114 
104.92852 
92.34394 
203.84649 
178.898 
44.20735 
246.80346 
103.79414 
0.27796 
31.05904 
264.72051 
336.72112 
280.01942 
199.22196 
326.65629 
276.38474 
185.01868 
105.58655 
107.63951 
21.96518 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
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2000622 
2000631 
2000633 
2000639 
2000642 
2000644 
2000651 
2000653 
2000658 
2000660 
2000661 
2000669 
2000673 
2000674 
2000675 
2000686 
2000695 
2000708 
2000714 
2000716 
2000720 
2000736 
2000737 
2000742 
2000749 
2000753 
2000770 
2000775 
2000782 
2000797 
2000800 
2000807 
2000811 
2000824 
2000839 
2000847 
2000851 
2000858 
2000864 
2000876 
2000883 
2000888 
2000897 
2000901 
2000925 
2000937 
2000939 
2000945 
2000951 
2000962 


2.416276 
2.7905149 
3.0197599 
3.0180846 
3.1960193 
2.6016046 
3.0236314 
3.0142155 
2.8545654 
2.5333023 
3.0162704 
3.0126289 
2.8147227 
2.9256113 
2.7677939 

2.588401 
2.5396531 
2.6712774 
2.5353112 
2.8129656 
2.8869458 
2.2019306 
2.5912677 
3.0109608 
2.2431611 
2.3290415 
2.2209539 

3.010628 
2.1799237 
2.5346964 
2.1930117 
3.0163759 
2.8951889 
2.7946416 
2.6138575 
2.7838528 
2.2282453 
2.8092647 
2.2082083 

3.011357 
2.2379184 
2.7086863 
2.5415771 
2.2237482 
2.6997365 

2.231828 
2.2465816 
2.6371401 
2.2093123 
2.9048251 


0.24156635 
0.0854026 
0.0877737 
0.1023474 

0.118776 
0.1544563 
0.0963724 

0.044892 
0.0617484 
0.1064984 
0.0366429 
0.0824955 
0.0107688 
0.1924971 
0.2042498 

0.269099 
0.1599762 
0.0832228 
0.0572934 
0.0860072 
0.0134885 
0.1649511 
0.2429092 

0.119093 
0.1736152 
0.2213364 

0.151461 
0.0748552 
0.0385567 
0.0602577 
0.2020272 
0.0668862 
0.0757905 
0.1328303 
0.1538848 
0.0939305 
0.0907124 
0.1035334 
0.1899383 
0.1077304 
0.1993118 
0.1942403 
0.0947401 
0.2215654 

0.081273206 
0.2178144 
0.1774328 

0.162079 

0.17412152 
0.1017686 


8.6414972 
18.93284 
10.90825 

8.57595 
8.14064 
1.04065 
10.76903 
11.28504 
1.50589 
15.21514 
9.25742 
10.78161 
2.87946 
13.51228 
9.80224 
15.68124 
13.85554 
3.48812 
14.27153 
8.4957 
2.3584 
4.37435 
12.36024 
11.21532 
5.38885 
10.08948 
4.3891 
9.28083 
5.26248 
4.50134 
4.2661 
11.30563 
3.13615 
8.11508 
12.60576 
2.48025 
2.39137 
8.88289 
5.44451 
11.3311 
4.71591 
13.85879 
14.32899 
3.44429 
21.068739 
3.69552 
2.5884 
32.84988 
4.1024509 
2.60193 


142.12698 
224.78721 
147.54125 
280.09687 
6.78999 
109.95042 
38.21433 
133.23306 
351.20658 
157.15003 
336.01605 
170.90917 
226.90118 
58.23889 
263.42953 
243.42817 
275.79136 
355.32832 
234.0746 
146.29921 
35.93197 
135.97225 
184.98656 
64.367 
109.88041 
61.48843 
44.8096 
298.00083 
80.53587 
238.47756 
325.26844 
132.34846 
130.95759 
141.74643 
338.27524 
271.19928 
141.24049 
67.30295 
163.21265 
151.14979 
285.70745 
124.23326 
258.05353 
265.30654 
299.73888 
243.80823 
327.28426 
318.39121 
253.21825 
145.66692 
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256.47956 
279.07028 
185.51247 
67.46954 
112.02326 
268.66057 
355.83247 
49.90575 
61.9558 
104.86147 
169.75515 
114.17126 
256.96914 
42.34252 
151.95144 
88.35887 
79.48882 
197.68247 
230.23879 
53.18106 
104.59844 
200.35955 
133.69582 
286.12007 
128.43565 
203.04489 
17.75077 
169.97933 
81.47485 
352.71322 
347.23211 
341.59535 
180.32103 
140.23895 
339.22259 
128.06198 
7.12275 
175.63874 
193.85566 
210.64758 
42.16645 
297 .67398 
22.97731 
68.06748 
201.94338 
72.00629 
5.94831 
161.12384 
129.49617 
225.13159 


228.26596 
175.92013 
261.16568 
300.6397 
51.50708 
230.99733 
296 .88522 
204.14329 
222.70176 
26.85243 
302.94086 
319.51092 
146.86493 
156.94179 
49.7925 
104.57222 
27.07343 
275.39522 
331.17093 
175.64261 
1.17256 
320.7405 
255.05591 
103.97612 
263.27234 
68.30372 
8.33105 
273.2585 
288.82101 
72.6863 
303.5713 
226.7289 
194.10173 
4745443 
144.00512 
166.21445 
115.0055 
73.08798 
326.6475 
313.11061 
232.0517 
349.84095 
293.86656 
201.12121 
171.13251 
311.97017 
215.58057 
7.3489 
113.11692 
69.99682 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
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2000963 
2000966 
2000968 
2000974 
2000975 
2001029 
2001043 
2001047 
2001052 
2001055 
2001058 
2001078 
2001079 
2001087 
2001088 
2001112 
2001129 
2001133 
2001140 
2001148 
2001185 
2001186 
2001215 
2001216 
2001223 
2001224 
2001245 
2001249 
2001252 
2001274 
2001286 
2001289 
2001306 
2001307 
2001314 
2001329 
2001336 
2001339 
2001350 
2001391 
2001401 
2001415 
2001416 
2001418 
2001422 
2001434 
2001442 
2001449 
2001500 
2001504 


2.2475028 
2.7208903 
2.8690829 
2.5327706 
2.8338758 
2.8902401 
3.0918753 
2.2407153 
2.2358075 
2.1983055 
2.1965147 
2.2700948 
2.8768912 
3.0136918 
2.2013686 
3.021494 1 
3.0275955 
2.1861828 
2.7718627 
3.0137128 
2.2379296 
3.0185508 
2.5776575 
2.2322234 

2.868483 
2.3050331 
2.8931802 
2.2241593 
2.6956455 
2.2289372 
3.0223851 
2.8612633 
3.1478297 
2.2510661 
2.2950475 
2.6180742 
2.8506578 

3.019179 

2.857687 

2.545694 
2.2263167 
2.2234386 
3.0236638 
2.2418907 
2.2478713 
3.0188203 
2.8732162 
2.2225604 
2.2420833 
2.3994278 


0.1378807 
0.1277792 
0.1349555 
0.111811 
0.0354774 
0.022333 
0.0468296 
0.1930689 
0.1440351 
0.2076132 
0.187668 
0.1382134 
0.0437337 
0.0953399 
0.1962884 
0.1014478 
0.0795809 
0.1868778 
0.1109592 
0.1155158 
0.1053667 
0.1082265 
0.1332289 
0.1793551 
0.0605204 
0.1986313 
0.077307 
0.0759426 
0.204297 
0.1135176 
0.0884699 
0.058457 1 
0.096918 
0.0962503 
0.1750754 
0.1709445 
0.0598015 
0.0560568 
0.0873096 
0.1678673 
0.1800958 
0.08697 
0.102367 
0.2038894 
0.1669438 
0.0609659 
0.0811964 
0.1423314 
0.190578 
0.1585168 


7.98932 
14.39219 
11.59647 

5.46304 

2.55965 

2.42891 

8.92769 

5.66428 

4.69483 

5.2722 

3.68964 

7.367 

1.17674 

10.0704 

7.65469 

8.99516 

8.60007 

5.37678 
14.13293 
10.84361 

5.70048 
10.75844 
15.91705 

7.60351 

2.55052 

7.87488 

2.88642 

4.87109 

33.8886 

4.39773 

9.73977 

1.61241 
14.91028 

3.94621 

5.24428 

14.4656 

3.19404 

8.68943 

2.93719 

7.58659 

7.28695 

3.42654 
10.04353 

7.19825 

2.67554 
10.81365 

1.25469 

6.63863 

7.44355 
11.04154 


62.56217 
72.62607 
209.00875 
86.76519 
38.85387 
30.1445 
159.56785 
78.33893 
99.67349 
147.21064 
221.93836 
93.94718 
329.6275 
30.51174 
54.58213 
302.99073 
269.61883 
58.33863 
72.1973 
145.68035 
71.99598 
43.21724 
123.81564 
121.67187 
41.07034 
258.27234 
151.89411 
259.11704 
141.08271 
327.29125 
200.93254 
193.25201 
274.51946 
233.97124 
264.77378 
132.21704 
97.51433 
291.08859 
139.67262 
103.57028 
277.67193 
329.4137 
353.02629 
355.20138 
201.72553 
152.79993 
221.18371 
110.83068 
20.04662 
94.96816 
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4.90994 
178.16174 
298.94054 
301.95786 

52.91719 
140.88708 
157.83728 
299.87517 
297 .39036 
176.30191 

93.99614 

43.88779 
106.23144 

28.59489 
319.43627 

86.60968 
134.63107 
306.58117 

311.2347 
174.52707 

2.23257 

294.1858 
265.72624 
144.58578 

10.09496 
128.83203 
207 .55894 
223.39854 

62.91168 
244 .48473 
102.70785 
117.98548 
134.41113 
207.16335 
144.00208 
165.03163 
218.75154 
166.96765 
239.65381 

85.14437 

70.88823 
240.54671 

62.48575 
324.07708 
170.75431 
142.70098 
126.91285 

131.9687 

16.77771 

51.38235 


44.92148 
144.49032 
78.47586 
106.25134 
345.36437 
327.05951 
246.86415 
162.4419 
80.72568 
0.29895 
318.1367 
68.71652 
163.95851 
336.51307 
74.07687 
259.12752 
180.36891 
287 .54725 
288.89636 
245.78469 
182.64599 
222.2892 
93.16039 
32.44269 
109.17591 
192.44571 
296.80457 
41.71452 
193.19854 
288.76466 
4.68589 
59.06832 
160.8342 
128.16734 
318.51489 
26.91872 
19.58499 
224 .50772 
327 .23494 
96.86838 
155.25928 
304.88173 
166.59297 
277.18144 
259.86018 
60.04093 
149.86304 
179.36625 
79.33813 
86.57975 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
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2001532 
2001533 
2001584 
2001601 
2001602 
2001619 
2001621 
2001636 
2001644 
2001648 
2001657 
2001665 
2001681 
2001707 
2001711 
2001717 
2001723 
2001755 
2001830 
2001842 
2001990 
2002000 
2002050 
2002052 
2002089 
2002090 
2002111 
2002156 
2002345 
2002411 
2002422 
2002430 
2002510 
2002830 
2012746 


Group 

2002062 
2002100 
2002340 
2003362 
2003554 
2003753 
2005381 
2005590 
2005604 
2033342 
2065679 
2066063 
2066146 


3.0037347 
3.0128199 
2.3765044 
2.2338815 
2.2446852 
2.2410031 
2.2300159 

2.234451 
2.5505794 
2.2356559 
2.3493066 
2.4134306 
2.6953244 
2.2187809 
3.0143118 
2.1959125 
3.0115989 
3.0902268 
2.1881582 
2.2662171 
2.1741611 
2.3822032 
2.3256838 
3.0071574 
2.5336949 
3.0715557 
3.0177442 
2.2423419 
3.0156037 
2.2254446 
2.3282179 
2.3627815 
2.2529475 
2.3780438 
2.2382462 


4 

0.9667013 
0.83206393 
0.84389537 
0.98946745 
0.97371428 
0.99774049 
0.94744835 
0.98568232 
0.92728217 
0.71851637 
0.91486343 
0.99093072 
0.78727129 


0.0553451 
0.0341753 
0.1942035 
0.12969 
0.1036564 
0.1755953 
0.1197028 
0.1280727 
0.153738 
0.2070075 
0.2345214 
0.2075508 
0.207426 
0.1709882 
0.1116368 
0.1286555 
0.0458738 
0.0498287 
0.0558756 
0.1801126 
0.0512001 
0.2970119 
0.2370149 
0.0836701 
0.156226 
0.1351973 
0.0909077 
0.2018135 
0.0784936 
0.0866982 
0.1986802 
0.2139148 
0.1963001 
0.2064344 
0.1908783 


0.18271178 

0.4364584 
0.44995917 
0.46856161 
0.28048062 
0.51478951 
0.29610804 
0.27958167 
0.40535123 
0.41798842 
0.26464041 
0.72009385 
0.48380682 


8.7858 
10.69007 
26.6425 
4.94351 
4.16426 
6.21376 
3.16935 
4.4334 
7.01216 
4.56623 
23.40366 
10.8324 
7.22492 
4.03832 
11.08177 
6.19124 
10.92085 
10.69528 
3.95267 
5.35433 
3.13114 
22.74792 
26.57894 
9.50768 
15.39516 
11.79893 
10.48847 
5.35454 
9.14384 
1.61434 
6.40677 
23.44678 
5.26923 
25.32361 
4.71175 


18.932519 
15.757434 
5.8539782 
9.9183703 
23.361458 
19.8093 
48.973045 
14.186004 
4.796818 
7.3416842 
1.2914011 
22.666222 
5.4099722 


330.92626 
156.87489 
305.46984 
74.77234 
75.1916 
61.59932 
182.00627 
168.50633 
270.95555 
130.50431 
105.44236 
91.66581 
94.62869 
6.28529 
134.93638 
340.64383 
150.00584 
157.33044 
147.57026 
153.5978 
193.75332 
292.18933 
72.61229 
213.98224 
102.80752 
340.05302 
167.35643 
17.28794 
304.07412 
131.05565 
160.1336 
46.00305 
102.98427 
49.12056 
265.31898 


108.62768 
170.87687 
211.51958 
152.50975 
358.67621 
126.29699 
58.562459 
216.34471 
312.00418 
82.004826 
178.35188 
351.92132 
102.31935 


82 


127.1195 
359.37335 
187.97207 
196.50973 

73.33204 
328.19457 
238.19504 
238.94113 
197.49709 
134.30702 

53.94026 

6.12477 
0.72186 

42.38917 
252.71894 
115.78748 

6.18978 
327 .98813 
335.13972 
125.43054 

12.25806 
129.81927 
170.76705 
207 .86004 
287 .08913 
337.43316 
233.75672 

4.19018 
139.06296 
129.58459 

52.03561 
309.62649 
209.36488 
140.44195 

81.63331 


147.94205 
355.99521 
39.938047 
54.982849 
359.38617 
43.743147 
37.413535 
34.416179 
82.422924 
167.28181 
14.927526 
151.06633 
84.562367 


284.23106 
324 .18222 
263.98458 
69.6855 
325.00527 
270.73248 
329.30148 
63.78235 
124.28199 
151.96224 
238.70086 
194.44845 
0.45966 
107.0253 
330.46795 
92.40988 
172.39803 
229 .63128 
22.26099 
264.22915 
241.19664 
140.16102 
36.08949 
301.1236 
59.54672 
329.75184 
284.67168 
170.84981 
287 .37178 
199.82879 
295.16215 
16.44569 
28.3665 
88.0689 
100.23076 


334.92171 

158.8253 
303.69494 
230.97821 

297.0219 
261.34727 
345.61516 
218.05496 
192.01877 
149.53901 
324.04271 

312.1291 
277.80779 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


WDWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW WW Ww 


AHRARHRHRHRPHRHRHRHRHAAA 


2066391 
2066400 
2068347 
2085770 
2085953 
2085989 
2086450 
2086667 
2087309 
2087684 
2088213 
2096590 
2099907 
2099942 
2105140 
3005972 
3005821 
3092114 
3012397 
3092124 
3005964 
3005969 
3005970 
3005973 
3092144 
3092156 
3009717 
3010201 
3010207 
3007848 
3011815 
3013030 
3013071 
3092192 
3081066 
3014184 
3021790 
3014113 
3014114 
3017039 
3015691 
3016523 
3027730 
3064315 
3017309 
3017060 
3020946 
3092226 
3019650 
3092245 


0.64231117 
0.85520974 
0.96282557 
0.99873332 
0.73876145 
0.88271351 

0.9680575 
0.85926966 
0.84745188 
0.85873344 

0.9539439 
0.90795755 
0.72853347 
0.92226308 

0.9142539 
0.89088169 
0.90800988 
0.93557259 
0.87624175 
0.68441915 
0.99299794 
0.75686007 
0.67082068 
0.90672544 
0.89755766 
0.97997322 

0.9132051 
0.93754452 

0.8655388 
0.83004916 
0.89669302 
0.87849063 
0.96311766 
0.73134122 
0.81647942 
0.70287493 
0.93250655 
0.81932815 

0.9032239 
0.87579594 
0.851 13036 
0.74123992 
0.87828555 
0.98185794 
0.91140605 
0.93722796 
0.81902244 
0.90718857 
0.67402486 
0.92884396 


0.68842863 
0.57256097 
0.37978055 
0.3449384 
0.70312537 
0.63293609 
0.41479377 
0.59473781 
0.46306937 
0.64269364 
0.59522617 
0.35118792 
0.59476 
0.1910585 
0.81702607 
0.66458077 
0.27176599 
0.09336139 
0.5514582 
0.50211708 
0.28400259 
0.39730327 
0.526364 1 
0.86843062 
0.28065043 
0.49926135 
0.36807412 
0.34638355 
0.20828137 
0.47478051 
0.3581831 
0.44056444 
0.31262107 
0.69855898 
0.64201038 
0.5042051 
0.34499429 
0.53002281 
0.50379849 
0.3181469 
0.44357995 
0.36744228 
0.73912509 
0.2096551 
0.11077323 
0.23454007 
0.46252854 
0.33269924 
0.66520869 
0.1120346 


38.891521 
9.0653163 
17.139554 
33.178463 
12.598508 
17.047146 
18.106551 
14.284786 
34.745734 
19.234294 
17.814259 
13.597697 
28.792699 
3.3313256 
32.514695 
7.2190179 
7.2538905 
12.377796 
2.0627017 
3.6341904 
23.752266 
6.8638137 
28.163343 
4.0352333 
3.8139909 
5.6585162 
31.677089 
12.773843 
16.718068 
25.492374 
16.208789 
3.3998126 
7.8072693 
30.350464 
53.296773 
2.9029512 
26.793225 
21.049912 
23.425461 
21.802114 
23.750628 
6.9691489 
13.432252 
7.2453491 
2.6226597 
6.5601162 
25.657332 
11.907748 
2.0164845 
10.776168 


244.9322 
79.918316 
245.97862 
18.403672 
180.54824 
130.29026 
124.91312 
208.40354 
294.31189 
162.11324 
114.31678 
75.804099 
225.62898 

204.46 
237.49863 
62.010882 

315.473 
329.18587 
2.7408046 
197.22215 
175.27593 
63.069077 
252.69674 
342.77703 
139.95032 
251.74134 
116.97656 
260.05979 
96.571071 
42.445929 
344.42506 
53.890616 
183.97257 
176.17726 
177.26469 
167.15519 
183.98793 
197.61161 
166.83655 
46.44937 
236.336 

74.57572 
280.10884 
85.995719 
313.34135 
327.40619 
155.93408 
240.07869 
80.792574 
232.31222 


83 


192.59839 
341.27615 

243.8006 
234.37669 
172.32944 
309.14389 
215.55604 
172.40027 
188.12379 
47.709401 
194.94151 
202.65036 
2.8083713 
126.39552 
281.51518 
193.51103 
336.42244 
354.30512 
253.64179 

179.0533 
349.68503 
205.68112 
356.52092 
322.39388 
150.27456 
55.816361 
141.55437 

203.7136 
16.625509 
180.85609 
356.81644 
309.00862 
260.72001 
359.82079 
359.48387 
6.2466136 
35.805081 
322.42211 
47.499929 
170.67926 
320.85172 
152.86115 

353.0528 
226.33812 
7.6536829 
299.48051 
253.33986 
223.55372 
9.8818082 

354.5409 


260.28645 
303.78632 
304.03906 
74.743373 
193.41695 
57.092802 
86.689225 
208.35792 
171.16429 
287.41396 
172.78666 
198.06344 
218.4728 
84.786507 
154.92805 
336.36256 
186.28381 
193.60748 
352.82705 
179.7846 
252.92589 
168.83864 
0.48117591 
99.988219 
206.58537 
72.906521 
290.00574 
229.81762 
314.20271 
86.169308 
222.69406 
327.02451 
53.391542 
118.70808 
110.45555 
188.83079 
63.591671 
146.69697 
92.326803 
55.367825 
257 .14348 
301.0636 
273.61263 
105.51781 
82.979128 
33.408243 
114.23065 
286.28021 
236.01441 
62.090196 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
49450 
54000 
49709 
54000 
54000 
50107 
50427 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
51077 
54000 
54000 
51081 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


ARR HRHRHRRHRPRPHRHRHRAPHRHRHRHRPHRHRPHRHRHRPHPHRHRHRPHRHRHRHRPHARHPHAPHRHRHRPHRAHRAPHRHRHRHRHRPARHRAAA 


3024030 
3029428 
3025763 
3025764 
3025765 
3031020 
3031703 
3028808 
3031176 
3031177 
3031178 
3092253 
3031183 
3031186 
3035165 
3035166 
3039898 
3036363 
3092260 
3042555 
3046648 
3053717 
3092272 
3061547 
3050515 
3054338 
3062815 
3054373 
3063789 
3063058 
3063823 
3057545 
3067492 
3067616 
3068066 
3071939 
3072413 
3069758 
3070801 
3072196 
3076722 
3072273 
3072291 
3076775 
3074756 
3079950 
3089251 
3092324 
3092325 
3092357 


0.89999766 
0.82941678 
0.85357291 
0.87831922 
0.74678009 
0.876503 
0.74122654 
0.86324958 
0.89556289 
0.83508053 
0.66184819 
0.81607747 
0.94696812 
0.93766812 
0.92871575 
0.92969018 
0.75784422 
0.81552756 
0.74399955 
0.91610235 
0.98475099 
0.87266158 
0.8258742 
0.95118128 
0.81138241 
0.92940488 
0.87028331 
0.88325596 
0.90366249 
0.87949364 
0.85448153 
0.91136333 
0.85678956 
0.8233799 
0.9402947 
0.8542365 
0.72538055 
0.71449449 
0.93967934 
0.98914235 
0.94421086 
0.98269502 
0.88860879 
0.87462622 
0.91403391 
0.93386971 
0.89155077 
0.95860272 
0.91409916 
0.90650678 


0.13956765 
0.55831156 
0.28634449 
0.41135064 
0.36001309 
0.89503754 
0.35849219 
0.42308796 
0.49535956 
0.56664333 
0.53291455 
0.46980368 
0.12351254 
0.32108762 
0.44688014 
0.45572011 
0.47657125 
0.43044014 
0.52412868 
0.11212565 
0.22113264 
0.37357638 
0.55130757 
0.31832272 
0.46688594 
0.16742376 

0.4223829 
0.24874524 
0.43874889 
0.26271471 
0.28864592 
0.78064631 
0.22462344 

0.3689377 
0.13738105 
0.17238766 
0.38255603 

0.4073083 
0.17625793 


0.057313022 


0.60426198 
0.02797089 
0.15914197 
0.49934979 
0.41209793 
0.32328704 
0.21371778 

0.240769 
0.34251754 
0.39299432 


1.6633228 
16.741231 
4.6958454 
2.6926995 
8.6005618 
25.676583 
6.8900577 
3.2721008 
11.560793 
13.775537 
26.470152 
3.8734528 
22.321499 
5.4192095 
40.260068 
28.578789 
32.143436 
2.6620264 
5.9877829 
2.8299469 
9.9847342 
51.180869 
19.643821 
16.093766 
10.355937 
0.89619274 
32.228534 
0.7818067 
11.744223 
17.412879 
7.6757121 
7.783026 
3.4909518 
17.513089 
5.7684736 
2.0261914 
8.1364752 
10.551364 
1.2921429 
11.62947 
20.663998 
5.2439975 
7.2897211 
23.74621 
5.209049 
12.090374 
50.207445 
13.877065 
7.5207799 
17.604042 


55.303545 
349.66292 
101.79881 
110.86572 
277.58834 
333.80753 
70.682966 
214.29152 
162.90292 
3.9767428 
155.80371 
345.26231 
331.20642 
178.31731 
358.52887 
52.953803 
358.1575 
55.249024 
30.318292 
72.769203 
304.64098 
124.4379 
177.92578 
312.84998 
350.69072 
14.741749 
29.823765 
237.8882 
33.888091 
50.862594 
55.931647 
69.374003 
97.371774 
297.87042 
115.62224 
122.57094 
109.55428 
331.00521 
31.947763 
357.82492 
22.148342 
183.09711 
189.24157 
32.651464 
205.39319 
295.94754 
317.70623 
170.77281 
154.15697 
11.3038 


84 


151.72637 
292.73797 
187.91717 
200.06223 
7.9287065 
324 .24248 
247.33611 
109.00818 
139.56794 
310.03182 
197.79925 
23.964009 
46.371289 
125.30451 
313.90075 
277.98821 
16.943112 
17.867756 
6.0866893 
7.7459151 
166.10892 
7.6049849 
354.37058 
211.71057 
224.30388 
131.29453 
187.37091 
293.04609 
228.71238 
229.98231 
221.95088 
213.58275 
189.32344 
30.585412 
242.83832 
195.57522 
234.13343 
353.49934 
342.489 
306.54524 
339.94382 
233.5423 
200.81708 
28.155709 
210.98323 
142.88185 
244.31983 
330.28055 
4.3166007 
213.96837 


221.34283 
117.89073 
3.8151398 
197.96191 
310.47211 

205.0929 
283.11028 
253.08418 
175.79773 
116.79868 
234.65129 

77.21737 
169.54532 
330.94723 
275.64076 
281.16831 
78.616265 
355.02601 
173.17855 
237.46661 
312.02601 

322.8322 
189.62089 
12.786194 
170.08819 
116.32414 
271.42008 
279.96198 
38.109895 

141.9183 
252.14788 
302.66343 
247 .32303 
332.21202 

185.9995 
217.07665 
173.96749 
284.26514 
166.51755 
94.254912 
205.64386 
352.59284 
356.38095 
337.96255 
203.67515 
155.09714 
122.67802 
356.41634 
117.71373 
38.677025 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
51624 
54000 
54000 
54000 
51662 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


ARR HRHRHRRHRPRPHRHRHRAPHRHRHRHRPHRHRPHRHRHRPHPHRHRHRPHRHRHRHRPHARHPHAPHRHRHRPHRAHRAPHRHRHRHRHRPARHRAAA 


3092339 
3092347 
3092370 
3092377 
3092380 
3092386 
3092390 
3102665 
3102680 
3102687 
3102718 
3102727 
3102728 
3102731 
3102744 
3102756 
3102762 
3102779 
3102787 
3114017 
3114023 
3114026 
3114104 
3117427 
3120884 
3117446 
3117447 
3117460 
3117468 
3120861 
3120863 
3125004 
3124996 
3125009 
3126183 
3127391 
3127401 
3127406 
3130459 
3131055 
3132092 
3133156 
3134264 
3134268 
3136734 
3137844 
3141527 
3141535 
3141538 
3143084 


0.94807271 
0.95414371 
0.91129913 

0.9619541 
0.87475176 
0.79674105 
0.88493377 
0.75103189 
0.79739141 
0.87176578 
0.67636878 

0.8799063 
0.77873106 
0.84083923 
0.85560809 
0.91555629 
0.99426257 
0.87507581 
0.97987696 
0.97888269 
0.86557884 

0.8198196 

0.8577063 
0.92124305 
0.92160627 
0.82348932 
0.76147828 
0.96703756 
0.98825788 
0.99226321 
0.76830954 
0.81884734 

0.7701479 
0.89868318 
0.92470305 
0.95500888 
0.71992412 
0.84475671 

0.9135579 
0.87651922 
0.93591413 
0.81697506 
0.96703422 
0.82550846 
0.90470877 
0.93000656 
0.79876331 
0.83739563 

0.7236569 
0.98456037 


0.49150795 
0.16607503 
0.525977 1 
0.48147436 
0.54598912 
0.77741935 
0.28666987 
0.3734296 
0.54626241 
0.38714096 
0.54206925 
0.54157069 
0.43775061 
0.38595466 
0.3736 
0.16340646 


0.013065778 


0.54652018 
0.1766903 
0.42840825 
0.22561468 
0.40303203 
0.36944345 
0.36305932 
0.56640928 
0.341939 
0.60181997 
0.30044155 
0.46273267 
0.52902156 
0.37557174 
0.39093655 
0.30533751 
0.26626348 
0.29777782 
0.37938048 
0.49553466 
0.31396358 
0.27422333 
0.43432919 
0.24289976 
0.69905755 
0.30996246 
0.28643033 
0.60068666 
0.12057551 
0.48723132 
0.43662198 
0.41051242 
0.20136283 


3.951064 
9.011133 
31.336012 
19.038812 
15.203483 
25.419665 
7.7052214 
18.179766 
27.159807 
30.995333 
4.8003637 
33.006701 
29.88594 
13.181771 
17.179764 
13.072181 
10.742631 
27.746329 
6.879208 
2.4597657 
3.1333566 
12.607732 
16.599775 
1.5475656 
30.14657 
6.5978745 
20.270204 
28.064829 
9.4897329 
22.554197 
10.628605 
40.85275 
4.3155863 
11.765305 
3.7633964 
2.9089622 
0.76027896 
6.1999594 
36.28443 
5.4175723 
6.905607 1 
12.748378 
13.538525 
1.3250455 
20.869143 
8.4784677 
30.700644 
9.6948341 
8.9801315 
2.331352 


327.3059 
13.218041 
27.399044 

12.90751 
57.856362 
30.332708 
25.592661 
239.71613 
69.743317 
92.982658 
306.08255 
294.73391 
287 .90059 
103.14846 
99.516019 
295.18464 
106.46935 

115.8071 
8.7510894 
81.442256 
137.63141 
137.70231 
234.34539 
347.22658 
262.94414 
164.19778 
203.66158 
7.9945393 
188.61692 
40.865148 
189.48132 
69.427745 
68.645093 
218.53494 
203.56944 
248.23483 
166.75667 
259.60686 
109.10504 
259.65972 
174.41506 
164.33324 
349.95955 
92.022439 
350.95397 
13.140614 
209.35617 
222.32902 
231.66157 
216.87263 


85 


262.35483 
241.34418 
208.76826 
199.00974 
135.96125 
212.38759 
133.00062 

358.3795 
208.48745 
219.67403 
318.32119 
25.324129 
323.84803 

149.5432 
205.13541 
306.60939 
100.60964 
147.21096 
331.56631 
272.77239 
210.35146 
217.99794 
94.099107 
57.667418 
25.307036 
223.15648 
148.24457 
31.452424 
226.69016 
8.3075997 
206.28661 
306.91181 
338.18528 
175.17419 
253.41942 
133.95352 
282.02629 
162.73931 
346.86651 
222.15883 
318.27888 
355.26495 
156.04858 
71.690414 
169.36365 
223.07557 
358.62917 
348.73063 
355.51291 
78.254887 


225.1088 
217.30782 
338.40154 
201.23957 
231.29666 
52.895742 
225.31406 
188.25934 
19.485949 
59.371016 
221.27778 
351.02314 
192.59727 

282.4749 
153.88991 

16.12534 
164.92244 
156.49332 
51.187861 
357.32802 
42.307944 
186.56678 
12.146551 
131.23247 
177.07735 
134.38613 
40.981494 
26.207003 
281.87596 

259.155 
4.8467488 
146.81126 
337.37454 
259.49608 
71.988181 
127.29616 

166.8299 
25.047051 
103.49846 
198.16622 
347.97534 
358.67083 
305.86572 
355.44719 
127.32394 
275.96035 
331.99547 
247.83159 

280.154 
72.980025 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
52242 
54000 
54000 
52271 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


ARR HRHRHRRHRHRPHRARHRAPHRHRHRPHRPHRHRPHRPHRHRAPHRHRHRPHRPRHRHRPHRHPHAPHRHRHRPHRARHRAHRHRHRHRHRAHRAAA 


3143121 
3144153 
3144155 
3144531 
3145517 
3146499 
3147579 
3150768 
3150774 
3151641 
3151644 
3151655 
3152309 
3152317 
3153508 
3153509 
3153530 
3154503 
3154513 
3154520 
3156302 
3157339 
3160723 
3160748 
3160799 
3160853 
3163736 
3164401 
3164404 
3164431 
3167348 
3167353 
3167367 
3170202 
3170242 
3170203 
3170204 
3170221 
3170208 
3172322 
3174187 
3175337 
3176187 
3177176 
3177188 
3177193 
3177197 
3177202 
3177226 
3177232 


0.90584107 
0.95516322 
0.90911752 
0.80953528 
0.87611211 
0.87482965 

0.9203212 
0.95752824 
0.93470009 
0.70741718 

0.8585099 
0.73081812 
0.89306537 
0.87232869 
0.84990274 
0.81317508 
0.81508903 
0.72744426 
0.94882172 
0.96054013 
0.85723098 
0.79330244 
0.88859382 
0.88518207 
0.82842354 
0.90779594 
0.77653201 
0.70092845 
0.78130319 
0.97915383 
0.90863991 
0.99058994 
0.88991985 
0.93049583 
0.89864059 
0.87876494 
0.93033425 
0.99414961 
0.96902587 
0.88403687 
0.96090346 
0.94402813 
0.88455581 
0.98950948 
0.88044837 
0.81459703 
0.90031102 
0.81796821 
0.88559948 
0.87856871 


0.23739246 
0.35943369 
0.21719253 
0.24208909 
0.38410512 
0.42612924 
0.11970839 


0.051791757 


0.24963774 
0.48594178 
0.39390777 
0.58168944 
0.21856995 
0.18199623 
0.380551 
0.27040152 
0.26213285 
0.5109053 
0.33060063 
0.14975606 
0.21046138 
0.49250313 
0.18325773 
0.30354972 
0.2097393 
0.31592396 
0.38182999 
0.59594346 
0.3368181 
0.14704531 
0.5445115 
0.12109868 
0.18053258 
0.19888123 
0.45058732 
0.26655962 
0.31303949 


0.013466245 


0.35503538 
0.22188425 
0.38494142 
0.40024487 
0.33025347 
0.27977013 
0.50520426 
0.26847244 
0.17108277 
0.28900036 
0.24938109 
0.49234168 


25.53355 
21.528325 
2.0854213 
34.160061 
7.3889096 
23.236807 
7.4793473 
15.269684 
13.220376 
23.375604 
13.050528 
6.6462784 
12.031238 
17.024551 
18.108305 
26.275789 
4.8810454 
23.505997 
23.901897 
10.795804 

0.63241404 
18.233456 
4.6515135 
3.6597329 
8.4615437 
25.432456 
12.146048 
16.508885 

3.794778 
17.058952 
28.539008 
2.5230027 
0.3560634 
19.460508 
29.268846 
5.7547748 
6.5240763 
4.3015637 

2.735198 
3.6171905 
17.667309 
23.023004 

5.141462 
4.6639542 
28.587726 
7.5899015 
7.9562364 


0.021373399 


6.7662088 
3.5104099 


245.93098 
265.95557 
160.69647 
81.700317 
301.85306 
286.83005 
139.95971 
346.05269 
178.04024 
177.59264 
21.684072 
359.37936 
196.34869 
199.60692 
70.456253 
30.512696 
39.08647 
215.53035 
250.048 
247 .33008 
215.81727 
124.69036 
178.77306 
176.57871 
274.20954 
200.72375 
220.16201 
22.958041 
188.87281 
212.81145 
57.597699 
42.306512 
65.466287 
89.811876 
86.451474 
89.91715 
99.18043 
265.1633 
86.544614 
299.09124 
300.91031 
157.35712 
336.708 
343.45054 
161.13348 
356.20573 
357.43335 
296.18068 
170.07059 
183.39904 


86 


351.2243 
317.58094 

118.5063 
178.61437 
23.557751 
43.944275 
172.95921 
23.978949 
167.88093 
196.51278 

339.2499 
29.591842 
181.81855 
168.68648 
306.68814 
19.241924 

351.8594 

203.9569 
154.67071 
238.12675 
210.40318 
311.55026 
316.20899 
47.968282 
326.31267 
355.12527 

321.8564 
200.59164 
59.273754 
287.43237 
140.64874 
225.08689 

138.0166 
165.60388 
138.05774 
222.81766 
134.48003 
83.895555 
127.98154 

28.23779 
45.800414 
216.07589 
50.031926 
55.822201 

202.3902 
337.56064 
343.24023 
31.319869 

196.2384 
142.00897 


346.24144 
263.35691 
274.21168 
246.44399 

298.6393 
283.29971 
223.42272 
75.424192 
155.74913 
158.54796 
342.06502 
322.20328 
202.52777 
275.52077 
322.20206 
24.021864 
56.376705 
308.37372 
73.518182 
294.20537 
248.98378 
87.736143 
88.522675 
321.24813 
130.70576 
314.04643 
329.76908 
116.86678 
206.13498 
269.39447 
280.45945 
97.336963 
5.1901772 
236.09578 
357.61805 
232.34244 
274.13873 
14.977607 
235.18049 
217.61998 
38.593947 
41.130168 
128.14063 
302.36837 
144.08901 
8.6811993 
181.87279 
14.154055 
197.01913 
280.59227 


54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
52724 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
52939 
52966 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


ARR HRHRHRRHRHRPHRHRHRAPHRHRHRHRPHRHRPHRPHRHRPHPHRHRHRPHRHRHRPHRPARHPHRPHRHRHRPHRARHRAPHRHPHRHRHRARHRAAA 


3177234 
3249978 
3179349 
3179363 
3180192 
3182186 
3182187 
3182823 
3182829 
3182833 
3183847 
3184475 
3249980 
3250193 
3250195 
3250293 
3251510 
3251512 
3252104 
3253645 
3254500 
3255174 
3255464 
3255465 
3255879 
3256321 
3256580 
3256583 
3257077 
3258062 
3258076 
3261401 
3261681 
3262331 
3262569 
3263232 
3263233 
3263448 
3263449 
3263451 
3263793 
3264188 
3264189 
3264547 
3265905 
3265909 
3266031 
3266035 
3267564 
3273458 


0.91350452 
0.82681962 
0.69653246 
0.78916702 
0.97997642 
0.95274176 
0.90121009 
0.83006288 
0.961 14304 
0.71192116 
0.91488208 
0.95062062 
0.95021121 
0.95199695 
0.94312385 
0.95092007 
0.90403261 
0.92044085 
0.96120257 
0.95427638 
0.87511021 
0.73778019 
0.86591757 
0.7670961 
0.90794321 
0.75098496 
0.89500976 
0.97640702 
0.94624489 
0.95420397 
0.9644237 
0.94051139 
0.94377277 
0.96831267 
0.85560093 
0.83752123 
0.88711985 
0.74868835 
0.76001344 
0.93149302 
0.64024467 
0.86840745 
0.84262456 
0.80883312 
0.872222 
0.88383727 
0.84659761 
0.94878176 
0.99121273 
0.89266375 


0.34920303 
0.39218262 
0.48849431 
0.59877773 
0.22337695 
0.56156984 
0.26577903 
0.40744127 
0.17028425 
0.49856275 

0.3513503 
0.56280186 

0.4167979 
0.17441825 
0.33803687 
0.12170564 
0.65680238 
0.35110701 
0.32880847 

0.1863123 
0.46481928 
0.41619723 
0.23763089 
0.42892764 
0.24989018 
0.44284208 

0.209597 1 
0.16192739 
0.25852582 
0.39678012 
0.22112866 

0.2439765 
0.16445192 

0.6554068 
0.16979483 

0.2980902 

0.1708655 
0.43382753 
0.40978133 
0.26650978 

0.7968061 
0.31327764 
0.24011097 

0.5361478 
0.56873666 
0.42113866 
0.30073285 
0.35599452 


0.068706638 


0.77031035 


33.467907 
4.1644361 
14.547286 
28.975304 
11.135007 
14.734804 
10.52491 
1.9081961 
35.070972 
22.055369 
25.453919 
29.334285 
14.339061 
56.22508 
6.2477966 
0.57411596 
15.91456 
5.9585478 
5.3345585 
22.053165 
21.334924 
18.418049 
18.701611 
4.7613942 
13.512538 
2.6217937 
19.353928 
36.517136 
14.059494 
3.7129817 
4.5080307 
16.218207 
1.2936976 
36.276971 
10.348597 
1.2031678 
12.187234 
3.0995667 
21.488824 
10.738706 
23.746543 
6.0667175 
12.086974 
28.149668 
5.8691155 
31.187501 
12.969629 
10.674391 
9.5265808 
16.320488 


195.47359 
191.24861 
115.75337 
203.02793 
214.71419 
235.26374 
101.99185 
243.48958 
78.132599 

79.22656 
82.670132 

263.9296 
327.20744 
151.86061 

75.46293 
175.15212 
119.11494 

173.8943 
199.45033 
356.7981 1 
46.654234 
180.36366 
302.13149 
202.37822 
14.576143 
48.487884 
155.44787 
12.606993 
203.10884 
29.862271 
211.98291 
225.83194 
233.53862 
208.92346 

54.94469 
285.35506 
253.39186 
307.15697 
85.680513 
120.95563 
122.70088 
263.47717 
265.61112 
269.94609 
117.66371 
116.00428 
296.89813 
113.41088 
147.02475 
344.39249 


87 


210.26718 
139.79199 
278.56009 
159.32865 
112.10254 
207.42161 
348.94624 
213.17629 
350.0398 
340.67122 
309.52078 
231.36764 
232.1766 
9.1122288 
109.19721 
28.533305 
65.154363 
333.85506 
280.81436 
226.8361 
94.378224 
359.49238 
233.48043 
322.48787 
200.48313 
136.67192 
77.702026 
138.64939 
56.491094 
120.62117 
294.21659 
297 .07694 
332.35792 
43.876742 
179.78693 
0.81535776 
15.930031 
302.8351 
157.51318 
115.66393 
130.7697 
47.285452 
7.8051078 
31.209738 
155.76728 
168.67007 
38.72732 
174.10611 
248.61695 
328.04398 


55.132809 
262.65627 
230.04016 
47.977419 

87.29048 
290.45178 
27.260244 
141.06958 
355.76643 
148.13048 
121.70682 
230.38682 
88.266623 
165.27237 
210.80311 
202.67682 
317.11312 

323.2308 
315.18661 
173.13211 
53.416337 
183.96342 
3.3642175 
207.26302 
239.92139 
208.67795 
230.72107 
238.04986 
130.79876 
303.80802 
275.85009 

286.7039 
220.29651 
42.361252 
299.55207 
244.86104 
209.35138 

107.2211 
99.054976 
228.79036 
23.609494 
167.70475 
264.56936 
243.68826 
265.87081 
191.20433 
170.98457 
131.31082 

326.2324 
297.39556 


54000 
53100 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
53241 
53244 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
53272 
54000 
54000 
54000 
54000 
54000 
53289 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
53433 


ARR HRHRHRHRHRPRPHRHRHRAPHRHRHRPHRRHRPHRHRHRAPHRHRHRPHRPRHRHRHRHPHRPHRHRHRPHRARHRPAPHRHPHAPHRHRHAHRAAA 


3273788 
3273782 
3274691 
3274905 
3275869 
3275978 
3276398 
3276601 
3276686 
3277400 
3278402 
3279867 
3283218 
3283227 
3283249 
3283679 
3283835 
3283950 
3285073 
3288855 
3288933 
3289173 
3289739 
3290865 
3291224 
3292261 
3293790 
3293831 
3293922 
3293923 
3297182 
3297356 
3297379 
3297628 
3297629 
3299721 
3304566 
3305028 
3306579 
3307228 
3307229 
3309039 
3309091 
3309092 
3309828 
3309832 
3309857 
3309858 
3311964 
3313739 


0.95945314 
0.76389989 
0.91854667 
0.93308734 

0.7532051 
0.77884535 
0.95886166 
0.79131389 
0.95138211 
0.85478599 
0.84051334 
0.98527456 
0.80376703 
0.85281863 
0.86343014 

0.7892811 
0.77941341 
0.86926394 
0.97601434 
0.89356017 
0.97558635 
0.99925501 
0.87749296 
0.98106367 

0.8655333 
0.84131167 
0.82682197 
0.94924874 
0.92386804 
0.83735026 
0.89288955 
0.93729026 
0.95814923 
0.74145111 
0.89095046 
0.94479231 

0.6716719 
0.95853738 
0.96807507 
0.84064591 
0.78419138 

0.7109288 
0.76032651 
0.81856961 
0.74387078 

0.8212713 
0.77172747 

0.9229554 
0.98261634 
0.91259519 


0.13537803 
0.38338784 
0.27307317 
0.33019919 
0.34020433 
0.38380306 
0.24585539 
0.39482129 
0.13570386 
0.33557573 
0.21493604 
0.79257712 
0.38177862 
0.29558723 
0.41102769 
0.49636097 
0.48307519 
0.42248981 
0.31997874 
0.36457579 
0.17580176 
0.30304443 
0.39370241 
0.28308734 
0.62371237 
0.41651256 
0.23366556 
0.37694757 
0.13387525 
0.22569769 
0.17018689 
0.56966238 
0.11600732 
0.42653472 
0.22521846 
0.23311402 
0.57517628 
0.15192594 
0.32044152 
0.26639019 
0.41420507 
0.55050252 
0.371895 
0.27264748 
0.33313323 
0.2963478 
0.32160915 
0.51225824 
0.64901708 
0.36592669 


30.002363 
20.484115 
12.946712 
3.7480425 
24.917136 
28.00374 
5.5684355 
9.1506832 
18.653347 
7.8992289 
2.9048383 
41.420794 
29.514513 
27.788514 
6.3146007 
10.639743 
6.0356351 
26.58707 
12.52351 
9.4605918 
3.9570556 
33.871509 
36.068592 
34.840215 
27.756769 
5.2041768 
25.622782 
5.005577 1 
2.4267858 
0.73313176 
5.6532049 
14.293632 
5.4166472 
24.545328 
0.2501009 
2.0866116 
23.030433 
21.623399 
32.651067 
17.248864 
16.853707 
19.581348 
12.795812 
3.6070987 
22.196585 
15.586801 
7.7314462 
14.137932 
4.1646477 
8.3160351 


329.84879 
353.00599 
0.23650582 
177.42193 
272.74031 
22.727001 
229.95645 
161.0855 
203.11495 
59.628165 
226.50305 
88.69818 
98.090608 
263.65459 
176.759 
289.8331 
114.63866 
120.89719 
309.81948 
48.218676 
334.9335 
155.09527 
165.88104 
23.566708 
354.90012 
8.67502 
14.945903 
195.25531 
346.0492 
196.83563 
32.652164 
58.989362 
216.10028 
223.74752 
37.454728 
49.351328 
69.435747 
288.1221 
75.70229 
84.890649 
282.3189 
288.69066 
274.41391 
71.203243 
282.7454 
281.79624 
100.74026 
127.80299 
123.37688 
305.98482 
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347.06403 
351.28138 
310.74842 
120.84034 
156.55353 
334.54877 
112.70238 
243.61206 
230.49457 
6.0811678 
181.14322 
42.781343 
349.71695 
190.52382 
250.09525 
194.62252 
0.58825253 
1.6808852 
204.75039 
108.6557 
119.6269 
54.330821 
320.37879 
205.51961 
228.05435 
151.88261 
185.79892 
304.934 
199.58911 
18.129628 
139.63883 
127.58469 
313.89196 
358.53892 
228.7358 
115.03663 
176.04117 
297 .85814 
238.56512 
149.93067 
9.2001035 
327.84503 
20.046772 
222.54974 
339.97538 
328.32981 
189.88574 
191.7352 
139.52651 
25.497647 


8.5373459 
192.40885 

150.4935 

160.1944 
293.86874 

268.0763 
75.266559 
140.13948 
316.55945 
42.263504 
102.53286 
197.64349 
85.132792 
9.5232994 
64.949008 
51.811201 
89.042008 
337.42373 
202.77323 
278.29264 
297 85286 
115.24323 
349.70934 
110.17821 
158.12306 
240.44876 
267.99709 
276.71907 
220.00037 
155.00452 
260.57564 
263.15198 
220.99061 
309.99123 
132.91306 
247.96876 
9.3177124 
119.57933 
87.195098 
220.28289 
165.71215 
281.67673 
181.31049 
148.62055 
249.41827 
221.01954 
184.97114 
14.612568 

173.2487 
43.203772 


54000 
53441 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
53647 
54000 
53652 
54000 
53655 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
53713 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 


ARR HRHRHRPRHRHRHRARHRAPHRHRHRHRPRHRPHRHRHRAPBRHRHRPHRHRHRHRPHRHRPHRPHRHRHRPHRHRHRAPHRHPHRHRHRAHRAAA 


3314789 
3315649 
3324253 
3324656 
3328632 
3329255 
3329278 
3330155 
3330538 
3330688 
3333079 
3337162 
3337325 
3338368 
3339082 
3341199 
3342642 
3342322 
3342323 
3343104 
3344169 
3347493 
3348144 
3350632 
3350633 


0.81990323 
0.67622812 
0.86368636 
0.69549487 
0.92250426 
0.95446956 
0.95908392 
0.84194719 
0.98525825 
0.84884937 
0.67190303 
0.83860414 
0.84773943 
0.80375651 
0.98674925 
0.95075741 
0.90638531 
0.80476877 
0.94928888 
0.88152606 
0.92191222 
0.98577785 
0.94018863 
0.94745779 
0.85107872 


0.4222185 
0.75501232 
0.33670256 
0.49762915 
0.34340855 
0.19860002 
0.17937767 
0.31558662 
0.26341245 
0.26080106 
0.58163186 
0.60517513 
0.57581932 
0.28451361 


0.046543451 


0.30076367 
0.23131827 
0.34642503 
0.28045436 
0.29156747 
0.32902583 
0.17475527 
0.39616796 
0.58305705 
0.46865622 


24.380862 
10.296226 
26.531002 
1.7812336 
14.622177 
1.5906721 
10.060842 
31.806967 
9.5375754 
16.418388 
5.9724883 

27.26727 

20.07954 

3.439796 
2.8266463 
1.4144836 

23.88335 
4.8343905 

5.865093 
4.5585453 
32.483792 
2.6419584 
11.569969 
5.5041004 
2.9133876 


146.84225 
303.38991 

162.3011 
336.43224 
15.228858 
280.91861 

183.9189 
35.992828 
232.82844 
34.794089 
41.039218 
291.85793 
115.27511 
4.9652624 

163.3305 

93.52479 

270.9998 

347.5677 
228.17418 
358.58023 
1.2876999 
194.39292 
195.46581 
225.42431 
92.154235 
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153.17094 
29.482862 
187.48359 
17.429135 
3.1455032 
154.81012 
242.88627 
317.68035 
212.56695 
22.616138 
17.670959 
174.50168 
29.325849 
124.76779 
332.95836 
110.27498 
261.22068 
182.91696 
305.53284 
165.97431 
224.38265 
38.266669 
315.13044 
299.84819 
68.436379 


103.18551 
217.22446 
57.179409 
140.24267 
356.30048 
276.65717 
281.19884 
246.80424 
246.14269 
331.26257 
14.966162 
318.38756 
199.74045 
235.02335 
225.00633 

136.8701 
144.12389 
190.23549 
190.56084 
205.79041 
103.38686 

110.3162 
265.13542 
272.59296 
235.83518 


54000 
54000 
54000 
53794 
54000 
54000 
54000 
53846 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54000 
54013 
54021 
54000 


ARR RHRARHRHRHRHRPHRHAPHAPHRHRHRHRHRHRPHRAHRAAA 
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APPENDIX B: 


Cost Function: 


function 


By Dr. 


Mh © 





MF © 


eventCost 
E10") ¥ 
runningCost 





primal 
primal] 


= primal 
primal 


Cost File. 
Pooya Sekhavat 
Last updated May 31 2007 by David Koeppel 





[eventCost, 


-nodes(1); 
.-nodes (end) ; 


-states (4,1); 
states (4,end); 


(t_£t - t_0)/m_f; 


runningCost] 


SPCost (pri 





0; 


Dynamics Function 


function residuals 


By Dr. 


























global CONST 

global EARTH 

global ast_ephem 

global SCALED 

% Spacecraft 

X_Sp = primal.states(1,:); 
y_sp = primal.states(2,:); 
Z_Sp = primal.states(3,:); 
m_sp = primal.states(4,:)j; 
vx_sp = primal.states(5,:) 
vy_sp = primal.states(6, 
vz_sp = primal.states(7, 
thrust_x = primal.controls 
thrust_y = primal.controls 
thrust_z = primal.controls 
thrust_mag = primal.contro 


r_mag 





Dynamics File 
Pooya Sekhavat 
Last updated May 31 2007 by David Koeppel 











primal.controls(5 


(1 
(2 
(3 
ls 

2); 


x 


x 


¥ 


x 
~~ wH 


OrbitDynamics (primal) 


x_spdot = 
y_spdot = 
z_spdot = 
m_spdot = 
vx_spdot = 
vy_spdot = 
vz_spdot = 

% Thrust 

6 Thrust 

% Thrust 
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mal) 





%-(m_f - m_0)/(t_f - 
primal.statedots(1,: 
primal.statedots(2,: 
primal.statedots (3,: 
primal.statedots (4,: 











statedots (5, 


Thrust magnitude 
Distance from sun 


l.statedots(6,: 
l.statedots(7,: 


M-FILES USED FOR FEASIBLE SOLUTION 


ls = primal.states; 




























































































.RU*%3)) .*x_sp./ 
































ED.RU*3)).*y_sp./ 

























































































ED.RU*3)).*z_sp./ 














| SP; 


U)).*thrust_y./m_ 


i) /SCALED.X) .*vx_sp; 
Ti) /SCALED.Y) .*vy_sp; 
i) /SCALED.Z) .*vz_sp; 


(r_m 





% Equations of Motion: 

residuals(1,:) = x_spdot - ((SCALED.Vx*SCALED.Ti 

residuals(2,:) = y_spdot - ((SCALED.Vy*SCALED. 

residuals(3,:) = z_spdot - ((SCALED.Vz*SCALED.Ti 

residuals(4,:) = m_spdot + 
(SCALED.TT*SCALED.Ti/SCALED.MU) .*thrust_mag./ (CONST.g*CONST.isp) ; 
residuals(5,:) = vx_spdot + 

( (SCALED .X*SCALED.Ti*CONST.mu) / (SCALED.Vx*SCALED 

3) - ((SCALED.Tx*SCALED.Ti) / (SCALED.Vx*SCALED.MU)) .*thrust_x./m 
residuals(6,:) = vy_spdot + 
((SCALED.Y*SCALED.Ti*CONST.mu) / (SCALED. Vy*SCAL 

3) - ((SCALED.Ty*SCALED.Ti) / (SCALED.Vy*SCALED.M 

residuals(7,:) = vz_spdot + 
((SCALED.Z*SCALED.Ti*CONST.mu) / (SCALED.Vz*SCAL 

3) - ((SCALED.Tz*SCALED.Ti) / (SCALED.Vz*SCALED. 
























































Events Function 



































function eventConditions = SPEvents (primal) 

& Events File. 

% By Dr. Pooya Sekhavat 

& Last updated May 31 2007 by David Koeppel 
global CONST 

global EARTH 

global ast_ephem 

global data 

global SCALED 

global AstrNo 

if (numel (ast_ephem) <1) 

ast_ephem = dlmread('ast_ephem.txt'); 

end 

% Spacecraft States x = [xX, y, Z, VX, Vy, vz] 
x0 = primal.states(1,1); x 
primal.states(1,end); 

yO = primal.states (2,1); Vy 
primal.states (2,end) ; 

z0 = primal.states (3,1); Zz 
primal.states (3,end) ; 

m0 = primal.states (4,1); 

vx0 = primal.states(5,1); vxf 
primal.states(5,end); 

vy0 = primal.states(6,1); vyft 
primal.states (6,end); 

vz0O = primal.states(7,1); vzft 
primal.states(7,end); 
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MU) ) .*thrust_z./m_ 





0 primal.nodes (1); 
f = primal.nodes (end); 


= primal.states (4,1); 
= primal.states (4,end); 





3 
rh © 
| 


oe 





= primal.parameters (1); 


0 Eccentric Anomaly 
f = primal.parameters (2); 





tO = primal.nodes(1)* SCALED.Ti+CONST.MinLaunchJD; 
tf = primal.nodes(end)* SCALED.Ti+CONST.MinLaunchJD; 




















M_Earth = EARTH.M0+sqrt (CONST.mu/EARTH.a%3) * (tO-EARTH.t0) ; SMean 
anomaly of Earth 

















theta = 2 * atan (tan(E0/2) / sqrt((1 - EARTH.e)/(1 + EARTH.e))); 
gamma = atan ((EARTH.e * sin (theta))/(1 + EARTH.e*cos (theta) )); 
SFlight path angle 
vr = EARTH.a * (1 - EARTH.e*%2) / ( 1 + EARTH.e*cos (theta) ); 
SMotion of Earth 

v = sqrt (CONST.mu * ((2/r) - (1/EARTH.a))); 

SMotion of Earth 






































SMotion of Earth converted into Cartesian Coordinates 
x = r*(cos(theta + EARTH.w) *cos(EARTH.lan) - cos(EARTH.inc)*sin(theta + 






















































































































































































EARTH.w) *sin(EARTH.1lan)); 

y = r*(cos(theta + EARTH.w) *sin(EARTH.1lan) + cos(EARTH.inc)*sin(theta + 
EARTH.w) *cos (HARTH.1lan)); 

Zz = r*(sin(theta + EARTH.w) *sin(EARTH.inc)); 

vx = v*(-sin(theta + EARTH.w - gamma) *cos(EARTH.lan) - 

cos (EARTH.inc) *cos (theta + EARTH.w - gamma) *sin(EARTH.1lan)); 

vy = v*(-sin(theta + EARTH.w - gamma) *sin(EARTH.lan) + 

cos (EARTH.inc) *cos (theta + EARTH.w - gamma) *cos (EARTH.1lan)); 

vz = v*(cos(theta + EARTH.w —- gamma) *sin(EARTH.inc)); 








X_init_ Earth = 
[x/SCALED.X, y/SCALED.Y, z/SCALED.Z,vx/SCALED.Vx, vy/SCALED.Vy, vz/SCALED.V 
Zee ce SInitial position and velocity of Earth 



































eventConditions = zeros(7, 1); 

eventConditions(1) = 1e5*(((x0 - 
X_init_Earth(1))* (SCALED.X/SCALED.RU))*2 + (yO - 
X_init _Earth(2))* (SCALED.Y/SCALED.RU))*2 + ((z0 
X_init_Earth(3))*(SCALED.Z/SCALED.RU) )%*2); 
eventConditions(2) = le5*(((vx0 - 


D.Vx/CONST.MaxVel 





X_init_Earth (4) ) * (SCALI ) 
X_init_Earth (5) )* (SCALED.Vy/CONST.MaxVel1) 
X_init _Earth(6))* (SCALED.Vz/CONST.MaxVel1) ) * 
eventConditions(3) = m0 - CONST.m0/SCALED 


























Hew od 

















ast_used = ast_ephem(906:910,:); 
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data =[]; 



























































data = CartesianAst (ast_used,Ef,tf); 

Sposvel = CartesianPV 
(tf,ast_used(:,2),ast_used(:,3),ast_used(:,4),ast_used(:,6),ast_used(:, 
5),ast_used(:,7),ast_used(:,8)); 

%[Delta_Vf,ind] = min(((vxf - data(:,4))* (SCALED.Vx/CONST.MaxVel)).*2 + 
((vy£ - data(:,5))* (SCALED.Vy/CONST.MaxVel)).*2 + ((vzf - 

data(:,6))* (SCALED.Vz/CONST.MaxVel)).%*2); 

[Delta_V,ind] = min(data(:,7)-X_init_Earth(7))j; 

eventConditions(4) = le-8*(((xf - data(ind,1))*SCALED.X)*2 + ((yf- 
data(ind,2))*SCALED.Y)*2 + ((zf - data(ind,3))*SCALED.Z)%*2); 
eventConditions(5) = 1le4*(((vxf - data(ind, 4))*SCALED.Vx)*2 + ((vyf - 
data(ind,5))*SCALED.Vy)*2 + ((vzf — data(ind,6))*SCALED.Vz)%*2); 
eventConditions(6) = M_Earth —- EO + EARTH.e * sin(EO); 
eventConditions(7) = data(ind,8) - Ef + data(ind,9) * sin(Ef); 





%SeventConditions (8) (Ef = 20) fmt 





AstrNo=ind; 


z KKKKKKKKKKKKKKKK 


END of FILE 








e5 


Path Function 


function PathFunction SPPath (primal) 


ol? 


Path File 
By Dr. Pooya Sekhavat 


oe 


ae 





















































global SCALED 

global CONST 

thrust_x = primal.controls(1,:);%.*SCALED. 
thrust_y = primal.controls(2,:);%.*SCALED. 
thrust_z = primal.controls(3,:);%.*SCALED. 
thrust_mag = primal.controls(4,:);%.*SCAL 
r_mag = primal.controls(5,:);%.*SCALED.RU 
X_Sp = primal.states(1,:);%.*SCALED.X; 
y_sp = primal.states(2,:);%.*SCALED.Y; 
Z_sp = primal.states(3,:);%.*SCALED.Z; 
PathFunction (1,:) = 

(thrust_y.* (SCALED.Ty/SCALED.TT)).*2 + 
(thrust_z.* (SCALED.Tz/SCAL 
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S$Gradual Cost Constraint 


KKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


Last updated May 31 2007 by David Koeppel 


























(thrust_x.* (SCALED.Tx/SCALED.TT)) .*2 








x & Thrust x direction 
vi % Thrust y direction 
2 & Thrust z direction 
Dg Pe %& Thrust magnitude 
, 





+ 


ED.TT)).*2 - thrust_mag.*%2; 


PathFunction (2,:) = (x_sp.* (SCALED.X/SCALED.RU)) .*2 
+ (y_sp.* (SCALED.Y/SCALED.RU)) .*2 + (z_sp.* (SCALED.Z/SCALED.RU)) .*2 = 
r_mag.%*2; 




















Main File 


% Spacecraft Initial Conditions in xyz 





oe 


Main File. 
By Dr. Pooya Sekhavat 
Last updated May 31 2007 by David Koeppel 


oe 


ole 


clear all; 


warning off MATLAB: fzero:UndeterminedSyntax 






































global CONST 

global EARTH 

global ast_ephem 

global data 

global SCALED 

global AstrNo 

CONST.AU = 1.49597870691e+08; % km SAstronomical Unit 

CONST.mu = 1.32712440018e11; % km3/s2 %Gravitational parameter of the 
Sun 

CONST.g = 9.80665e-3; %Skm/s2 SStandard acceleration due to 
gravity 

CONST.isp = 4000; oS sspecific impulse 

CONST.vinf = 3.5; % km/s shyperbolic excess velocity 
CONST.DayToSec = 86400; % s SOne day in seconds 
CONST.MinLaunchJD = 57023*CONST.DayToSec; SS SProblem 
start date 00:00 01 January 2015, (MJD) 

CONST.MaxArrivalJD = 64693*CONST.DayToSec; % s sProblem 
end date 24:00 31 December 2035, (MJD) 

CONST. year = 365.25 * CONST.DayToSec; SS SOne year 
in seconds 

CONST.DtoR = pi/180; smultiplicative conversion 
from Degrees to Radians 

CONST.DisTol = 1000; % km 

CONST.vTol = le-3; % km/s 

CONST.m0O = 1500; & kg sinitial mass 1.5 ton 
= 1500 kg 

CONST.mProp = 1000; % kg Sprop mass 1 ton=1000 
kg 

CONST.Tmax = 0.le-3; % kg.km/s2 Smax thrust 











ton.AU/Day*2= 0.1(N=kg.m/s%*2) 


ole 














SCALED.RU = 10*CONST.AU ; % km = 15*AU (Since a's are given in AU 
untes) 3 

SCALED.X = 10*CONST.AU ; 

SCALED .Y 1O*CONST.AU ; 
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SCALED.Z = 10*CONST.AU ; 

SCALED.Ti = 4*sqrt (SCALED.RU*3/CONST.mu) ; 

SCALED.Vx = 4*sqrt (CONST.mu/SCALED.RU) ; 

SCALED.Vy = 4*sqrt (CONST.mu/SCALED.RU) ; 

SCALED.Vz = 4*sqrt (CONST.mu/SCALED.RU) ; 

SCALED.MU = 1000; 

SCALED.Tx = (CONST.mu*SCALED.MU) / (SCALED.RU%2) ; 

SCALED.Ty = ee RU%2); 

SCALED.Tz = (CONST.mu*SCALED.MU) / (SCALED .RU%2) ; 

SCALED.TT = pee enon MU) / (SCALED .RU%2) ; 

%1* (SCALED . MU*CONST.isp*CONST.g) / (SCALED.Ti) ; 

EARTH.a = 0.999988049532578*CONST.AU; % (km) 

EARTH.e = 1.671681163160e-02; SEccentricity 
EARTH.inc = 0.8854353079654e-03 * CONST.DtoR; SInclination in 
RADIANS!!! 

EARTH. lan = 175.40647696473 * CONST.DtoR; SLongitude of 
Ascending Node in Radians!!! 

EARTH.w = 287.61577546182 * CONST.DtoR; sarg of perihelion in 
RADIANS!!! 

EARTH.MO = 257.60683707535* CONST.DtoR; SMean Anomaly at 
epoch in RADIANS!!! 

EARTH.t0O = 54000*CONST.DayToSec; SEarth epoch time 
(MJD) 

CONST.MaxPos = 15*CONST.AU; %5.909335501968675e+008; S15*CONST.AU; &% 
km 

CONST.MinPos = 0.1*CONST.AU; 

CONST.MaxVel = 100; % km/s 

Xmin = 

CartesianPV (CONST.MinLaunchJD, EARTH.a, EARTH.e, EARTH. inc, EARTH.w, EARTH.1 
an, EARTH.MO, EARTH.t0); 

6 % Xmax= 

CartesianPV (CONST.MaxArrivalJD, EARTH.a, EARTH.e, EARTH.inc, EARTH.w, EARTH. 
lan, EARTH.MO,EARTH.t0O); 

% Define the problem: 

JPL.cost = 'SPCost'; 

JPL.dynamics = 'SPDynamics'; 

JPL.events = 'SPEvents'; 

JPL.path = 'SPPath'; 








Setup the knots and some more constants 








8*CONST. year/SCALED.Ti; 
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(CONST.MinLaunchJD-CONST.MinLaunchJD) /SCAL 





ED.Ti; 


tfiMax = (CONST.MaxArrivalJD-CONST.MinLaunchJD) /SCALED.Ti; 
tfMax_low = 8*CONST.year/SCALED.Ti; 























tfGuess = tfMax_low+ (t f{Max-tfMax_low) /2; 

knots.locations = [tOmin tfGuess]; % time free problem 
knots.definitions = {'hard', 'hard'}; 

knots.bounds.lower = [tOmin tfMax_low]; & [tOmin 
tfiMax_low]; 

knots.bounds.upper = [tOmin_up tfMax]; % upper bound was tf0 = 
fixed. 

knots.numNodes = [15]; 

% Bounds 





oe 


bounds.lower.states 
-CONST.MaxPos/SCA 
CONST .MaxVel/SCAL 
CONST .MaxVel/SCALED.Vz]; 
bounds.upper.states = [ CONST.MaxPos/SCALED.X; CONST.MaxPos/SCALED.Y; 
CONST .MaxPos/SCALED.2; CONST.m0/SCALED.MU; CONST.MaxVel/SCALED.Vx; 
CONST .MaxVel/SCALED.Vy; CONST.MaxVel/SCALED.Vz]; 


[-CONST.MaxPos/SCALED.X; -CONST.MaxPos/SCALED.Y; 
ED.Z; (CONST.m0-CONST.mProp) /SCALED.MU; —- 
D.Vx; -CONST.MaxVel/SCALED.Vy; —- 

















EF] el 





[t] 



































6) 
a 
i 


bounds. lower.contr [-CONST.Tmax/SCALED.Tx; -CONST.Tmax/SCALED.Ty; 
-CONST.Tmax/SCALED.Tz; 0; CONST.MinPos/SCALED.RU]; 
ia 












































bounds.upper.controls = [CONST.Tmax/SCALED.Tx; CONST.Tmax/SCALED.Ty; 
CONST.Tmax/SCALED.Tz; CONST.Tmax/SCALED.TT; CONST.MaxPos/SCALED.RU]; 
bounds.lower.events = [0; 0; O; O; O; O; O]; 

bounds.upper.events = [0; le5*CONST.vinf*2/CONST.MaxVel*2; 0; le- 





8*CONST.DisTol*%2; le4*CONST.vTol*%2; 0; O; 0.09]; 


bounds.lower.path = [0; 0]; % convexifying the nonconvex cone; 
lower bound = - inf 

bounds.upper.path = [0; 0]; 

bounds.lower.parameters = [0; 0]; 





bounds.upper.parameters [40*2*pi; 50%*2*pil; 


oe 





ae 


Provide a guess 


oe 





























guess.states(1,: = [Xmin(1)/SCALED.X, CONST.MaxPos/SCALED.X]; 
guess.states(2,: = [Xmin(2)/SCALED.Y, CONST.MaxPos/SCALED.Y]; 
guess.states (3,: = [Xmin(3)/SCALED.Z, CONST.MaxPos/SCALED.2]; 


= [CONST.m0/SCAL 
; 
Xmin (4) /SCALED.Vx, CONST.MaxVel/SCALED.Vx]; 
Xmin (5) /SCALE CONST .MaxVel/SCALED.Vy]; 
= [Xmin(6)/SCALED.Vz, CONST.MaxVel/SCALED.Vz]; 


guess.states (4, 


) 

) 

) 

:) D.MU, (CONST.m0- 
CONST .mProp) /SCALI 

:) 

2) 

=) 








al 
1w) 
ls 
Gq 


guess.states (5 
guess.states (6 
guess.states (7 



































Hee 
< 
Kg 
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guess.control 


s , CONST.Tmax/SCALED.Tx]; 
guess.controls 
s 


0 
0, CONST.Tmax/SCALED.Ty]; 
guess.control 0, CONST.Tmax/SCALED.Tz]; 
guess.controls(4,:) = [0, 
norm([CONST.Tmax; CONST. Tmax; CONST. Tmax],2)/SCALED.TT]; 
guess.controls(5,:) = [norm([Xmin(1) ;Xmin(2);Xmin(3)],2)/SCALED.RU, 
norm([CONST.MaxPos; CONST.MaxPos; CONST.MaxPos],2)/SCALED.RU]; 











Myc) 
2,:) = 
372) 








[ 
[ 
[ 
[ 














guess.time = [tOmin tfGuess]; 

guess.parameters (1) = [15*2*pi]; 

guess.parameters (2) = [20*2*pi]; 

beginCPU = cputime; 

[cost, primal, dual] = dido(JPL, knots, bounds, guess); 
RunTime = (cputime - beginCPU) /60 


cost 


Feasibility (primal); 


oe 


Bootstrap 


oe 


ole 


guess.states = primal.states; 

guess.controls = primal.controls; 

guess.time = primal.nodes; 

knots.locations = [primal.nodes(1), primal.nodes (end) ]; 
knots.numNodes = [60]; 


% 


JP AP WP ol? 


ole 





ole 


beginCPU = cputime; 
[cost, primal, dual] = dido(JPL, knots, bounds, guess); 
RunTime = (cputime - beginCPU) /60 


ol? 





ae ole 


ole 


Cost 


ole 


ol? 


Feasibility (primal); 
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